Interaction corrections to the thermal transport coefficients in disordered metals: 

quantum kinetic equation approach. 
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We consider the singular electron-electron interaction corrections to the transport coefficients in 
disordered metals to test the validity of the Wiedemann-Franz law. We develop a local, quantum 
kinetic equation approach in which the charge and energy conservation laws are explicitly obeyed. 
To obtain the local description, we introduce bosonic distribution functions for the neutral, low- 
energy collective modes (electron-hole pairs). The resulting system of kinetic equations enables us 
to distinguish between the different physical processes involved in the charge and energy transport: 
elastic electron scattering affects both, while the inelastic processes influence only the latter. More- 
over, the neutral bosons, though incapable of transporting charge, contribute significantly to the 
energy transport. In our approach we calculate on equal footing the electrical and thermal conduc- 
tivities and the specific heat in any dimension. We found that the Wiedemann-Franz law is always 
violated by the interaction corrections; the violation is larger for one- and two-dimensional systems 
in the diffusive regime Tt <C fi and it is due to the energy transported by the neutral bosons. For 
two-dimensional systems in the quasi-ballistic regime Tt ^ fi the inelastic scattering of the electron 
on the bosons also contributes to the violation. 

PACS numbers: 71.10.Ay, 72.10.Bg, 72.15.Eb 



I. INTRODUCTION 

It is well known that the measurement of the thermal 
transport coefficient may provide additional information 
about the scattering processes in disordered metals. In 
particular, the Wiedemann-Franz^ law holds as long as 
elastic scattering"^^ dominates in the system: 



(1.1) 



where k, and a are respectively the thermal and electric 
conductivities in the system, T is the temperature in en- 
ergy units (fcs = 1) and e is the electron charge. On the 
other hand for the deeply inelastic forward scattering the 
Wiedemann-Franz law is violated^, so that the Lorentz 
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number L is smaller than the universal value, L < 

Historically, the transport (in particular thermal trans- 
port) coefficients were ffist being calculated using the 
Boltzmann equation (BE)^. The advantage of this ap- 
proach is that it allows for the clear separation of the 
scale in the problem: a particle moves freely most of the 
time and rarely scatters on other particles or impurities. 
The BE is applicable on the time scale much larger than 
the time it takes for the scattering to happen, so all the 
scattering events are encoded into the local collision inte- 
gral. All the quantum mechanical part of the calculation 
is reduced then into the solution of the scattering prob- 
lems for the relevant physical processes. This gives the 
precise form of the collision integral but it does not affect 
the general structure of the BE. The great advantage of 
the BE is that its structure illuminates the relevant con- 
servation laws. 

In the late 50s an alternative approach became popu- 
lar - the so called Kubo formulas'*. In this approach the 
transport equation is not derived but rather the connec- 



tion of the transport coefficient to the equilibrium corre- 
lation function of certain current operators is used. (The 
Kubo approach to the thermal transport was claimed to 
be put on rigorous footing by Luttinger^ based on the as- 
sumption that there exists some spatial scale in the sys- 
tem so that for perturbations smooth on that scale the 
gradient expansion is possible). Being exact, the Kubo 
formulas are formally applicable even in the regime where 
the transport equation can not be justified (separation of 
the evolution into free motion and rare collisions is not 
possible). 

However, in practice, the possibility of explicit calcula- 
tion within the Kubo formula is somewhat limited. The 
most spectacular results of the Kubo-formula calcula- 
tions - such as Maki-Thompson®'^, Aslamazov-Larkin^ 
and weak localization^ corrections to the electrical con- 
ductivity - require a small parameter which is the same 
parameter that determines the applicability of the Boltz- 
mann equation. It means that all those effects can be 
also described in terms of quantum corrections to the 
collision integral (for weak localization it was done in 
Ref. 13). The most relevant for this paper effect - the 
Altshuler-Aronov^'^ interaction correction to the electri- 
cal conductivity in two dimensions^*'"'^'* 



5 (J A A = 



ln( A 



(1.2) 

originates from the elastic scattering of the electrons on 
the self-consistent potential (Friedel oscillation)*^'^^, and 
it can be once again obtained from the correction to the 
collision integral*^. 

The success of the Kubo formulas in the description of 
the quantum and interaction effects in thermal transport 
is by far more modest and controversial. Particularly, 
despite a 20-years history there is no consensus on the 
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answer to a natural question: how does the logarithmic 
correction to the conductivity (1.2) translate into a cor- 
rection to the Wiedemann- Franz law (1.1)? 

A first attempt at answering this question was made 
by Castellani at al.'^^ by analyzing Ward identities for a 
disordered Fermi liquid; they found that the Wiedemann- 
Franz law should hold for interacting disordered elec- 
trons. Their claim was later disputed by Livanov ct 
al.^'': in a "quantum kinetic equation" approach,"^® a 
logarithmic divergence for the thermal conductivity in 
two dimensions was found to have even the sign contrary 
to the Wiedcmann-Franz law. More recently Nivcn and 
Smith^^ applied Kubo formula and found again a loga- 
rithmically divergent contribution (for the Coulomb but 
not the short-range interaction) in addition to the one 
that follows from the Wiedemann- Franz law. 

The reason for this confusion in the literature is 
twofold. Technically, the identification of the correct 
form of the current operator is complicated by the pres- 
ence of the electron-electron interaction (the energy cur- 
rent operator in the form defined by Luttinger^ is cum- 
bersome to use due to the presence of the additional dis- 
order and intc;raction potentials in it, whereas the su- 
perficially more elegant expression in the Matsubara fre- 
quency representation in fact does not correspond to any 
conservation law for the interacting system and violates 
gauge invariance - see Appendix B). Physically, the use 
of the diagrammatic calculation within the Kubo formula 
prevents one from a clear identification of the relevant 
scattering processes, because each diagram taken sepa- 
rately describes some mixture of such processes and does 
not have its own physical meaning. 

This situation calls for the development of the kinetic 
equation description, which takes into account the inter- 
action correction of the Altshuler-Aronov type both for 
the electric and thermal transport. The advantage of this 
approach is that it makes possible to keep track of the 
conservation laws explicitly and thus excludes any am- 
biguity in the definition of the currents. This paper is 
devoted to the development and application of this ma- 
chinery. 

We will use units with h = 1 throughout the paper 
and we will restore the Planck constant in the final re- 
sults only. The remainder of the paper is organized as 
follows: in Sec. II we discuss some general features of the 
kinetic equation approach using a simple "toy model". 
In Sec. Ill, we present our final expression for the ki- 
netic equation describing interacting electrons in disor- 
dered metals. Section IV summarizes the results for the 
thermal conductivity and the specific heat obtained by 
solving the kinetic equation. The derivation of the kinetic 
equation is presented in Sec. V, while the calculation of 
the transport coeflicients and the specific heat is given 
in Sec. VI. Some mathematical details are relegated into 
Appendices. 



II. STRUCTURE OF THE KINETIC 
EQUATION: CURRENTS AND SPECIFIC HEAT 

The purpose of this section is to show how the structure 
of the kinetic equation permits the proper identification 
of the relevant currents. We will remind how to calculate 
the specific heat from the kinetic equation once the con- 
servation laws are obtained (this enables a direct check 
against the much simpler thermodynamic calculation). 
We will discuss the locality requirement for a proper ki- 
netic equation. The latter requirement will dictate the 
number of necessary degrees of freedom (i.e. independent 
distribution functions) which should be introduced into 
the kinetic description. 



A. Kinetic equation and conservation laws 

As a concrete example, we consider here electron-like 
and hole-like excitations coupled to neutral bosons in the 
presence of an external electric field E. (As we will see 
later, a system of interacting electron can be eff'ectively 
described at low temperatures by such a coupled sys- 
tem for the scattering at small momentum transfer in 
a particle- hole channel.) The kinetic equations for elec- 
trons and bosons have the form: 



— + vpn ■ V -I- evpn ■ E— 
at oe 



f = Ste{f,N} (2.1a) 



dt 



■ v{Lv)n ■ V 



N = Sh{f,N} 



(2.1b) 



where / = /(e, n; t, r) is the distribution function for the 
electrons with charge e, vp is the Fermi velocity, and n is 
the direction of the momentum. Energy e is counted from 
the Fermi level so that /(e > 0) describes the electron- 
like excitations and 1 — /(— e), e > corresponds to the 
hole-like excitations. Having in mind only singular in T 
corrections, we neglect the dependence of the electron 
velocity on energy (electron-hole asymmetry) as it does 
not introduce anything but small correction regular in 
powers of T^. 

The bosonic function N = N{u, n; t, r) is the distribu- 
tion function for the bosons with velocity v{uj). All the 
interaction effects are included into the collision integrals 
Ste and Stb; for example, an electron-like excitation can 
decay into a less energetic electron and a neutral boson, 
or an electron and a hole can annihilate into a neutral 
bosons, etc. By locality, the collision integrals depend 
on the same variables as the distribution functions, i.e. 
Ste = Ste{e, n; t, r) and Stb = Sti,{uj, n; t, r). 

In thermodynamic equilibrium and E = the Fermi 
function for fermions and the Planck function for the 
neutral bosons 



fF{e) = 



exp(e/T) + l' 



Np{u)) = 



exp(a;/T) - 1' 



(2.2) 
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solve the kinetic equation. The temperature T here is 
a constant determined by the initial conditions for the 
kinetic equation. 

Being the effective description for the slow dynamics 
of the original quantum mechanical system, the kinetic 
equation must respect the conservation laws of the orig- 
inal system: (i) total charge conservation and (ii) total 
energy conservation. Those two conditions are enforced 
by the requirements for the collision integrals 



and 



/ 



dev{SU{f,N})r^ = Q, 



(2.3a) 



and 



j de e ,y{Ste{f,N})r, + J doj u b{u){Stb{f,N}}r. = 0. 

(2.3b) 

Here is the density of states (DoS) of the electrons 
(we will neglect its energy dependence) and b{Lu) is the 
density of states of the bosons. We also introduced the 
short-hand notation for the angular integral 



(2.4) 



where $1^ is the total solid angle in d dimensions. 
Let the electron density p be given by: 

p{t,r) = ep J de{f{e,n;t,r))n. (2.5) 

Integrating Eq. (2.1a) over the energy and the direction 
of the momentum and using Eq. (2.3a) we arrive at the 
continuity equation 

dp 



dt 



+ V • j = 0, 



(2.6) 



with the electron current density defined as 

j{t,r) = eyvF j de{nf{e,n;t,r))n- (2.7) 

(Strictly speaking Eq. (2.6) fixes only the longitudinal 
component of the current - an arbitrary curl may be 
added to Eq. (2.7). We will not consider the effect of the 
magnetic field here, so we will disregard such magnetiza- 
tion currents.) 

We turn now to the analysis of the energy conservation. 
We multiply Eq. (2.1a) by ve and integrate over e and n. 
Next, we multiply Eq. (2.1b) by Lob{Lu) and integrate over 
the U! and n. Adding the two results together we obtain, 
with the help of Eqs. (2.3b) and (2.7) 

+ V ■ 7^ 

' ^ Jtot 



dt 



j-E 



(2.8) 



where 



(2.9a) 
(2.9b) 



utot =Ue{t,r) +Ubit,r)\ 
Ue{t,r) = u j dee {f{e,n;t,r))n; 

Ub{t,r) = J duj ijj b{ijj){N{ijj,n;t,r))n, (2.9c) 



3 tot ~ 3e "I" 3bi 



(2.10a) 

jl{t,r) =vfv j dee (n/(e, n; r))„; (2.10b) 
jl{t,r) = j duujb{u)v{u){nN{LO,n;t,r))n(2.10c) 



The right hand side of Eq. (2.8) is nothing but the 
Joule heat. For the homogeneous system the gradient 
term in the left hand side vanishes, and by virtue of the 
energy conservation the expression (2.9) must be iden- 
tified with the total energy density of the system. On 
the other hand, ior E = Eq. (2.8) has the form of 
the continuity equation for the energy density; therefore 
Eqs. (2.10) has to be identified with the total energy cur- 
rent density. This statement is not entirely trivial. One 
could imagine that for an interacting system the DoS 
entering into the expressions for the charge (2.5) and for 
the energy density (2.9) are renormalized differently. The 
energy conservation equation (2.8) excludes such a pos- 
sibility. 

The conservation of energy (2.8) is valid for any rate 
of the energy flow in and out of the system. On the other 
hand, the collision integrals in Eqs. (2.1) define a certain 
time scale Tin - the dynamics slow on the scale of Tin 
can be characterized by the distribution functions (2.2) 
with the temperature dependent on time, T{t) [correc- 
tions to such an adiabatic description are of the order 
of Tindt InT]. Substituting such form of the distribution 
function in Eqs. (2.9), and the result in Eq. (2.8), we find 
for the homogeneous system 



dT 

cAT)^=3-E, 



where 



d_ 



f de efpie) + dio u b{u!)Np{u!) 



(2.11) 



(2.12) 



is nothing but the specific heat of the system. The latter 
quantity may be calculated independently by applying 
the standard diagrammatic technique for the equilibrium 
systems. The agreement of such a calculation with the 
structure of the kinetic equation result (2.12) will be the 
most important check of the consistency of our descrip- 
tion of the thermal transport. 



B. Locality of the kinetic equation and the number 
of degrees of freedom 

The local in space and time form of the collision in- 
tegrals is clearly a simplified description. Actually the 
collision integral may be nonlocal on the time scale of 
the order of h/T and on the spatial scale of the order 
of Tivp/T. We will call such a description local, and the 
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description where the non-locahty is involved on larger 
spatial and time scales - non-local. 

The number of distribution functions to be introduced 
into the description is governed by the locality of the ki- 
netic equation. Let us use the model of Eqs. (2.1) to 
illustrate the point. We had a local description in terms 
of the fermionic and bosonic distribution functions. One 
can, however, try to eliminate the bosonic distribution 
function and obtain a description in terms of the elec- 
tronic degrees of freedom only. 

Assuming that the deviation of the distribution func- 
tion from its equilibrium value is small, one can linearize 
the bosonic collision integral to the form 

SU{f,N} = -i(N-N{f}) (2.13) 

where I is some positive definite integral operator and 
N{f} is the functional of the fermionic distribution func- 
tion /(e), such that for /(e) = /^(e), N{lo) = Np{u). 
Using Eq. (2.13) one can formally solve Eq. (2.1b) 



N : 



IN{f} 



(2.14) 



Substituting Eq. (2.14) into Eq. (2.1a) we apparently ob- 
tain the kinetic equation in terms of the electron distri- 
bution function only 



— -I- vpn ■ V -I- evpn ■ E— 
at oe 

1 



5t:{/} = su f, 




(2.15) 



If we are interested in the linear response to a weak and 
smooth external perturbation, the description in terms 
of this single kinetic equation is completely equivalent 
to the original coupled system (2.1). However, there are 
clearly drawbacks: the presence of the integral operator 
/ in the collision integral makes it non-local on the scale 
determined by the kinetic equation itself rather than by 
the temperature. Moreover, though it is still easy to 
identify the continuity equation for the electron charge 
using Eq. (2.3a), there is no longer a relation similar to 
Eq. (2.3b). This is why the analysis of the energy con- 
servation law becomes cumbersome: the calculation of 
specific heat and energy current requires the time expan- 
sion of the collision integral, which in turns seems to re- 
quire the knowledge of the concrete form of the inelastic 
collision integral. 

The example we have just considered is somewhat 
trivial because the separation of the system into the 
fermionic and bosonic modes was given from the very be- 
ginning. The problem we are dealing with in this paper 
is how to include the collective modes of the interacting 
electron system into the kinetic equation. Indeed, in this 
case any calculation will give the result in terms of the 
electronic distribution function only, and it is not clear a 



priori how to introduce the occupation numbers for the 
collective modes into the description. 

As we will show, it may be possible to reverse our 
previous argument. We will consider a system of inter- 
acting electrons and we will find that the interactions are 
described by a non-local collision integral. Therefore we 
will introduce bosonic degrees of freedom that will enable 
us to rewrite the non-local kinetic equation in terms of 
coupled, local kinetic equations. Then we will be able to 
identify the energy density and energy current density as 
sums of fermionic and bosonic contributions. The con- 
crete example will be discussed briefly in the following 
subsection. 



C. Degrees of freedom for the kinetics of 
disordered Fermi liquid 

We now focus on the disordered, interacting Fermi liq- 
uid. For simplicity, we consider the interaction in the 
singlet channel only. Our goal is to show that the ther- 
modynamic result for the interaction correction to the 
specific heat has indeed the kinetic equation structure 
(2.12). As a result, we will be able to determine the nec- 
essary number of bosonic degrees of freedom for the local 
kinetic equation. For the paper to be self-contained, we 
briefiy review the thermodynamic approach, referring the 
reader to the literature^^ for further details. 

The thermodynamic calculation of the specific heat Cy 
is based on the relation between Cv and the thermody- 
namic potential O: 



QT2 



(2.16) 



The thermodynamic potential can be written as the 
sum of the thermodynamic potential for non-interacting 
quasiparticle Jlo and a correction SVl associated with the 
soft modes in the system. Keeping such a correction is 
legitimate because it turns out to be a more singular 
function of temperature than the correction due to 
the electron-hole asymmetry. 

The correction 5^1 \s given by the sum of the so-called 
ring diagrams, see Fig. 1. The Matsubara representation 
for this diagram is 



59. 



(217) 



Here F is the interaction constant, a;„ = 27rTn are the 
bosonic Matsubara frequencies and 11 is the polarization 
operator. The explicit expression for this operator is not 
important for the present discussion and will be given 
later, see Eq. (5.43a). 

A straightforward calculation, relegated to Appendix 
A, enables us to rewrite Eq. (2.17) as 



5Vl . 



dbJ f 1 , w \ 
— - coth — 
27r \2 2TJ 



{2-kY 



ImTr 



Inr^-ln^s 
(2.18) 
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the collective modes to the energy density in the non- 
equilibrium case 

POO 

Ub= duJoj[NP{uj)b''{uj)-NS{uj)¥>{oj)], (2.20) 

^0 

where A'^'' = = Np in the equilibrium and they have 
to be found from some kinetic equation otherwise [this 
definition requires that Eq. (2.8) holds for an arbitrary 
distribution frmction]. A similar expression can be ob- 
tained for the contribution due to the interaction in the 
triplet channel by introducing the additional propagator 
and distribution frmction N'^ . This means that the 
proper local kinetic equation must include four distri- 
bution functions: one for fermions, /(e), and three for 
bosons, NP''^'3{lo). In the subsequent sections we will be 
able to derive such a description. 



FIG. 1: Leading singular contribution to the thermodynamic 
potential for the clean system. The shaded box corresponds to 
F/v, defined through the two particle vertex , see Ref. 26; 
the solid lines are coherent parts of the electron Green's func- 
tions. For the disordered system, the polarization bubbles 
should be dressed by the impurity scattering^^. 



The explicit expressions for the bosonic propagators C, 
CP are not relevant [they can be found in Eq. (3.11) - 
the trace should be understood as a sum or integration 
over all variables other than tj, g]; we just mention here 
that = in the absence of interaction, F = 0. Sub- 
stituting Eq. (2.18) into Eq. (2.16) and performing an 
integration over w by parts, we find 

d f°° 

Scv = ^ dLJUjNp{(j)[¥{(j) -b'^{oj)], (2.19a) 
c*-* Jo 

where the densities of states are defined as 

bP{u;) = -Im / ^5^Trln£^; (2.19b) 

TT J (27r)" 

b3{u) = -Im. f -^d^TvlnCB. (2.19c) 
TT J (27r)" 

The function b''{uj) has the physical meaning of the den- 
sity of states (DoS) of the bosonic degrees of freedom 
in the system (soft electron- hole pairs). The function 
b^{co) has the meaning of the density of states of ficti- 
tious bosons (we will call them "ghosts" ) which describe 
the soft electron-hole pairs in the absence of interaction. 
The physical meaning of the minus sign in front of h^(uj) 
is that with the formation of the collective modes some 
degrees of freedom are removed from the description of 
the non-interacting system; the introduction of the ghost 
bosons in the last term in Eq. (2.19a) takes this reduction 
into account. 

Comparison of Eq. (2.19a) with Eqs. (2.12) and (2.9) 
hints the following expression for the contribution of 



III. FINAL FORM OF THE KINETIC 
EQUATION AND SCATTERING PROCESSES 

This section summarizes the final form of the quan- 
tum kinetic equation, the conservation laws and the cor- 
responding currents. The derivation of these results is 
presented in Sec. V. 

In accord with the previous section, the kinetics of 
the system is described by the electronic distribution 
function f{e,n;t,r), the "distribution functions" of the 
bosonic singlet and triplet excitations, NP and N'^, and 
the "distribution function" of the ghost excitation, . 

The electron distribution function f{e,n;t,r) is diag- 
onal in the space of the momentum directions. On the 
contrary, the bosonic excitations are characterized by the 
density matrices N°'{uj, q; rii, nj;t, r) [a = p, a, g] which 
may be not diagonal in the space of the momentum di- 
rection n. Only in the thermal equilibrium 

fegie,n;t,r) ^ fpie) , 

N"g{u;, q; n„ nf.t, r) = 0^5 {nlrij) Np{lu) , 

with the Fermi and Plank distribution functions given by 
Eq. (2.2), the matrices N°'{uu,q]ni,nj;t,r) acquire the 
diagonal form'"*^. However even out of equilibrium these 
matrices have the property 

7V"(a;,gr;ni,nj) 

= - [TV" (-W, -q; nj, Ui) -I- r2<j(5 (n^j)] . 

(hereafter, the spectator t, r variables might be sup- 
pressed.) 

Strictly speaking, f{e,n;t,r) is the 2x2 density matrix 
in the spin space and iV^ is the 3x3 density matrix in the 
angular momentum L = 1 space; however this will not be 
important in further calculations and we will write the 
equations for the diagonal components only. To account 
for the three folded degeneracy of the triplet mode we 
will explicitly introduce factors of 3 in the corresponding 
collision integrals and currents. 
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For the sake of compactness, we will use an operator 
notation for matrices in the space of momentum direc- 
tion, so that for example N should be understood as an 
operator acting on a function a{ni) as follows: 



Na 



dvt " 

ini)= I -^N{n,,nj)a{nj). (3.3) 



The kinetic equation for the electrons in the electric 
field E (we will not consider the magnetic field effects) 
has the canonical form: 



dt + v-V + ev-E 



f{e,n;t,r) = Ste(e, n; i, r); 



Ste = St,/ + St""!/, NP} + 3SC " {/, N-} 
-4St"'{/,7VS} + St"'{/}, 



(3.4) 



where the first term on the right hand side is the "bare" 
collision integral 



Str{ni,nj) = 



6 (n^n 



(3.5) 



with 6ij = Uirij, and the other terms, which will be 

written shortly, take into account the interaction effects. 

The bosonic distributions, for a = p,a,g, are governed 
by: 



^^-^■,dtN"^ + [s"iiO,qy,VN"} (3.6) 



+i 



h: 



);iV-] 



St {N°',f}{uj,q;ni,nj;t,r) 



where the commutator and anticommutator are defined 

as 

^A;B^ = ^{AB + BA); A;B = AB - BA. (3.7) 



The operators ff^^ acting in the angular (momentum 
direction) space are defined as 



HZhi^^Q) = vq 



and the velocity operator is 



1 + F" 



(3.8) 



= + . (3.9) 



dq 



dq \1 + F° 



The action of the operators in the angular space is the 
same as in Eq. (3.3) and they have the following forms: 

F» = 

[F-]{n,,nj) = F^{e,,) (3.10) 
[F'']{ni,nj)=uV{q)+FP{0ij), 



where FP'''{0) are the Landau Fermi- liquid interaction 
parameters. The angular independent vV{q) term takes 
into account the long range part of the Coulomb density- 
density interaction. 

To characterize the density of states for the bosonic ex- 
citations we introduce the propagators jC"{u),q;ni,nj), 
a = p,a,g as 



C" = 1, 



(3.11) 



They describe the propagation of the electron-hole pair 
which is scattered by the disorder potential. This prop- 
agation is affected by the corresponding interactions for 
a = p,a, whereas it reduces to the usual diffusion for the 
ghosts. 

We are now prepared to write down the conservation 
laws which must be satisfied by the collision integrals in- 
dependently of their explicit form or the particular shape 
of the distribution functions. The conservation of the 
number of particles is guarded by the condition 

J St {f,N°'}{e,n;t,r)dnd€ = 0. a = g,p,a; 
j St (e, n; t, r) dn de = 0, (3.12a) 

and the impurity collision integral (3.5) conserves the 
number of particles on each energy shell 



/ 



Stx/(e, n; t, r) dn = 0. 



(3.12b) 



The conservation of energy during purely electron-elec- 
tron collisions is dictated by 



/ 



eSt {f} {e,n;t,r) dnde = 0. 



(3.12c) 



Finally, the conservation of energy during the electron- 
boson collision is guaranteed by the conditions 



+ /TV 



eSt" {/,7V«}(e,n;t,r) 



dn de 



£"(a;)St " {/,7V"}(a;;t,r) 



dtu 
2^ 



Tr 



4«hM;>C"(a;) 7V«(a;;t,r) 



(3.12d) 



jduj 



2tt 



for a = g,p, a, where the trace is defined as 



TtAB 



f dnidn2 f d'^q 

J ^2 J j^A{q;ni,n2)B{q;n2,ni) 



(3.13) 

The existence of the conservation laws (3.12) imme- 
diately enables us to establish the expressions for the 
conserved currents in the spirit of Sec. 11 A. We find by 
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integrating both sides of Eq. (3.4) over e and n: 



p{t,r) =ei' j f{^,n;t,r) 



(3.14) 



de dn 



j{t,r) = evvF j nf{e,n;t,r) 



dedn 



which express the conservation of charge in terms for the 
usual charge density and electric current density - cf. 
Eqs. (2.5)-(2.7). 

Turning to the energy conservation, we multiply 
Eq. (3.4) by e and then integrate over n, e. Similarly, 
we pro-multiply Eq. (3.6) by take the trace (3.13) 
and integrate over u). Adding the results together, one 
finds: 



ot 



dt ' ' 

Utot = Ue+Up + 3Mct - 4?ig 



(3.15a) 



The electronic contributions to the energy density and 
current density are given by 

Ue{t,r) = v j ^^^^ef{e,n;t,r) 

, . f de dn \ /„ , n 

3eit^r) = i^VF J -^enf{e,n;t,r). (3.15b) 

The contributions of the bosonic neutral excitations are 

Ua{t,r) = J Tr| ^^^^^ £"(a;)jV"(a;;t,r)^ '^'^ 

3l{t,r)= j Tr{s"(w)£"(w)iV"(a;;t,r)} 



27r 

u] du) 



2tt 
(3.15c) 

for a = g, p, a. 

Equations (3.14)-(3.15) constitute our main results: 
the conserved currents are defined in terms of the distri- 
bution functions of the quasiparticlcs which describe the 
low energy excitations of the interacting electron gas for 
interaction in the particlc-hole channel. In contrast with 
previous calculations^""^'', wc explicitly show the valid- 
ity of the continuity equation for the energy transport; 
no such proof has been presented before in the quantum 
kinetic equation approach'^*. Moreover we believe that 
the form of the energy current in those references is not 
correct, since it is not gauge invariant - see Appendix B 
for more details. As an additional benefit, our approach 
enables us to clearly identify the contributions of the col- 
lective modes and the scattering processes involved this 
last task is accomplished by analysing the explicit form 
of the collision integrals, which is is also needed to cal- 
culate the transport coefficients. The detailed derivation 
of the collision integrals is presented in Sec. V; here we 



summarize the results and give them physical interpreta- 
tion. 

To shorten the formulas we introduce the following 
combinations of the distribution functions: 

"^trM {^,'^,Q;t,r) = iV" {u;,q;ni,nj;t,r) 

x[f{e,nk;t,r) - /(e - w, nfe; r)| (3.16a) 

+QdS{rMnj)^^f{e,ni;t,r) [1 - /(e - a;,nfc;t,r)] | 
which, as it follows from Eq. (3.2), has the property 

jderfj.kl{^,u;,q) = j deT%.ik{e,-uj-q) (3.16b) 
and 

^ij-M (e, ei;w) 

= /(e - rii) [1 - /(e, Uj)] /(ei, nfe) [1 - /(d - m)] 

- f{e,ni) [1 - f{e-uj,nj)]f{ei-uj,nk) [1 - /(ei,n;)]. 

(3.16c) 

It is easy to check that T = \1/ = in the thermal equi- 
librium (3.1). The last combination enters the collision 
integral in the symmetric form 



1 r 

4 



(3.16d) 



It is worth noticing that the terms involving four distribu- 
tion functions / are in fact cancelled from ^f*; moreover 
it has the properties: 



J (e,ei;a;) rfedei =0 
a(w) / e'if'lj {e,ei;uj) dedcidco = 



(3.16e) 



for any even function a (a;). Finally we introduce the 
vertex 7 for the impurity scattering 



1 



-fid [6 {rijUk) - 5 {niUk)] 



(3.17) 



T(%) 

and the shorthand notation 

The explicit expression for the boson-electron collision 
integral is then 



St {L0,q;ni,n2;t,r) 



X (luTf 32;4i (e, t^. 9; t, r) + ^li'^is-M (e- 9; r") } , 

(3.18) 

for a = g, p. a. The formula for the electron-boson col- 
lision integral can be conveniently decomposed into local 
(Z) and non-local (n) [in the sense of Sec. II B] parts 



St 



St," 



St. 



(3.19a) 
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a) 



h h 



e-h 





FIG. 2: The scattering amplitudes leading to the creation of 
the same electron and hole out of different electron hole pairs 
(double lines) (a-b) and their interference contribution (c). 
The impurity is denoted by filled circle. 



The local part of the coUision integral is 

(3.19b) 



where the bar indicates Hermitian conjugation 



(3.19c) 



Using Eq. (3.11) and the definitions (3.13), (3.17) we can 
verify that the pair Eqs. (3.18)-(3.19b) satisfies the en- 
ergy conservation law Eq. (3.12d) on its own. 

Equation (3.19b) satisfies the particle number conser- 
vation law (3.12a) as well. To check this, we make the 
change of variables (w, q) {—(^, —q) in the terms con- 
taining and then use Eqs. (3.16b) and (3.19c) to 
rewrite the integral of Eq. (3.19b) in terms of only: 



St; (e, ni)dedni = 



drii . . . dn^^ 



X {7?2[>C^4T4"i;2i(e,c.,g) - /:?4T4"3^i2(e,c.,q)] 

+ 7I1 [£?4T?2;2l(e, ^, Q) " ^24T4%;12(e, ^, O)] } 

We perform the integration using the delta functions 
in Eq. (3.17) and we obtain the result antisymmetric with 
respect to the ni <-!■ n2 permutation. Hence the above 
expression vanishes after the ni^2 integrations. 

The physical meaning of the collision integrals (3.18) 
and (3.19b) is the following. In the absence of the dis- 
order the electron-hole pair propagates for an infinitely 
long time. Due to the impurity potential the decay of the 
pair into the electron and hole moving in different direc- 
tions as shown in Fig. 2 is allowed. Equations (3.18) and 
(3.19b) are the probabilities for such a decay. [See also 
Sec. VD after Eq. (5.67) for further discussion.] 



The non-local contribution to the collision integral 
2 f du! I f d'^q f dn2 . . . dn-j 



St„ (e,ni) 



duj 1 
2^w2 



V 

X7?37l6['C?4-'Ci4]sinh 



(27r) 

2 [ ^ 
2de 



[fie, Tie) - /(e,n4)] 
X y dei{ [T^7.g4(ei,cj,g) + T^7.46(ei,tJ,q)] £?2 

+a7[Tf?2;13(ei,^,'7)+T?2;3l(ei:'^,g)] } 

(3.19d) 

by construction obeys its own conservation law 

e™St^'"(e,ni)de = 0, to = 0,1, (3.19e) 



i.e. it conserves the energy and the number of the elec- 
trons moving for a given momentum direction n. More- 
over, one can see that the collision integral (3.19d) does 
not contribute to the linear response at all because in 
the thermodynamic equilibrium T = and / does not 
depend on the angle. The non-locality of this collision 
integral indicates that the task formulated in Sec. II B 
has not been quite accomplished. Technically, this non- 
locality can be decoupled by introducing a density matrix 
non-diagonal in the boson-ghost space. We choose not to 
pursue this line because the term (3.19d) does not con- 
tribute to any observable quantity we are interested in 
and does not affect any conservation laws. 

The electron-electron collision integral can be split into 
elastic, non-local and local parts: 

St (e,ni) = Stgj (e,ni) -|-St„ (e,ni)-l-St; (e,ni). 

(3.20a) 

The elastic term describes the scattering of the electron 
on the static self-consistent potential created by all the 
other electrons 



~e-e 2^ f dio 1 
St., = -Re / 



d'^q f dn 



. dri: 



625 
-713746 



X [CP + 3£^ - AC%^ [/(e - ng) - /(e - m)] 



'Cf4/(e,n3) + ClJ{t,ni) 



(3.20b) 



Its physical origin is discussed in details in Ref. 17. Being 
elastic, it conserves the number of particle for each energy 
shell: 

St^,''(e,ni)dni = 0, (3.20c) 

as one can see from the property — —7!]^ of the vertex 
(3.17). 

The nonlocal term 



St, 



-Re 



dw 1 
2^w2 



d'^q 
{2tj:Y 



dn-). 



. dn 



dei{[C 



£§4] [LP + 3£" 



(3.20d) 
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describes the inelastic electron-electron collisions during 
which bosons and ghosts act as virtual states. [The func- 
tion were introduced in Eq. (3.16d)]. The real part 
being an even function, we can use Eq. (3.16e) to verify 
that Eq. (3.20d) obeys the conservation law 



/ 



e'"St„";(e,ni)de = 0, m = 0, 1. 



(3.20e) 



As indicated, the same law is satisfied by the local (and 
elastic) term: 

r> e-e 2 fdto 1 f cfiq f dn2 ■ ■ ■ driQ /" , n 

= u J 2^Zr^ J (2^^ J j '^'^^^^^^ 

X [f(ei,n,)[l- f{ei-co,n3)] + {n,^n3)]. (3.20f) 

Therefore equation (3.20e) enables us to conclude that 
both collision integrals (3.20d) and (3.20f) do not affect 
the transport coefficients [when they can be considered 
as perturbations in comparison to the bare impurity col- 
lision integral]. 

We note that, although it might not be evident, the 
present form of the kinetic equation permits the proper 
identification of the inelastic kernel that determines the 
phase relaxation time - further details can be found in 
Appendix G. 



IV. SUMMARY OF THE RESULTS FOR 
THERMAL TRANSPORT AND SPECIFIC HEAT 

In this section we present our final answers for the in- 
teraction corrections to the thermal conductivity and the 
specific heat. They are obtained by solving the kinetic 
equations and then substituting the solutions into the 
definitions (3.15) for the energy and energy current den- 
sities. The explicit calculations are performed in Sec. VI. 
We will consider the short range impurities t{6) = r. We 
report our results for quasi one- dimensional and three- 
dimensional systems in the diffusive limit Tr <C U; for 
two-dimensional systems, we will not put such a restric- 
tion on the temperature range. However common to all 
dimensionality is the zeroth harmonic approximation for 
the Fermi liquid constants [see Eq. (6.14)]. 



A. Thermal conductivity 



In the absence of the magnetic field the thermal con- 
ductivity tensor is diagonal, = KSfj,^, and we will 
write the expression for the diagonal components only 



The first term is given by the Wiedemann-Franz law 
Kwj? = LaT with the inclusion of the interaction correc- 
tions to the conductivity and the Lorentz number given 
in Eq. (1.1). The second term causes a violation of the 
Wiedemann-Franz law. In the diffusive limit and for low 
dimensionality the main contribution to An is due to the; 
long range nature of the bosonic energy transport, which 
originates from the long range part of the interaction in 
the singlet channel. In the quasi ballistic case a large 
contribution comes from the inelastic scattering of the 
electron on the bosons as well. Smaller corrections arise 
due to the triplet channel bosonic transport and to the 
energy dependence of the elastic scattering. 

For quasi one-dimensional and three-dimensional sys- 
tems in the diffusive limit we write 

Ak = 5kP + iSn" + 5Kei 

where the bosonic corrections include the ghost contri- 
butions 



K", a = p,a 



[see Eq. (6.12) for the definition of k"] and we neglect 
the inelastic contributions 6Kin which are smaller by the 
parameter Tr/h. 

For quasi one- dimensional systems the explicit expres- 
sions are: 



X 1-1+3 



^ /3 



DT 



(4.2a) 

(4.2b) 



DT ,1 fnoe 

——akhi^ 
h \ T 



/27r 



DT 



VT+Ff-1 



(4.2c) 



where a is a length of the order of the wire width and 
k = v47re^ is the inverse screening length in the bulk. 
For three-dimensional systems the results are: 



5Ke 



(4.3a) 



5kP = 




as: 



In these expressions, ((x) is the Riemann zeta function 
and C(3/2) « 2.612 , C(5/2) « 1.341. 

For two-dimensional systems, we separate the correc- 
tions due to the singlet and the triplet channel interac- 
tions: 



Ak. 



(4.1) 



Ak = Aks + 3AKt 



(4.4a) 
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The singlet channel contribution is, with logarithmic ac- 
curacy: 

Aks = —01 2n— In — ^ 



^92 (n- j In (^1 + ^ 



n'TfTT\\ (Ef 



(4.4b) 



where k = 27re^i^ is the two-dimensional inverse screening 
length. The cross-over functions gi and §2 are given in 
Eq. (6.42). Here we note that <7i(a;), .92(2;) ~ 1 for x ^1 
and gi{x) ~ S/x, 52(2;) — 14x^/15 for x » 1. 
For the triplet channel we have: 



T 
■ 18ft 



l-^ln(l + Fo-) 



Tr <C1 



Tr > 1 

(4.4c) 

In the diffusive limit Tt <^ h our results are consistent 

with those of Ref. 22, even though the form of the energy 
current operator used in this reference is, to our opinion, 
incorrect - see Appendix B. 



B. Specific heat 



The specific heat is given by 



Cv = —yT + 5Cv 



(4.5) 



where the first term is the usual non-interacting elec- 
tronic contribution and the second term is the bosonic 
interaction correction one. 

For quasi one- dimensional and three-dimensional sys- 
tems 

d 

5cv = (1 + 3Ad) ' aa (4.6) 

The two terms in the first bracket are respectively the sin- 
glet and triplet channel contributions. The singlet chan- 
nel term is considered in the unitary limit and is therefore 
independent of any interaction parameter. On the other 
hand, the Fermi liquid parameter for the interaction in 
the triplet channel enter into Eq. (4.6) as 



Arf = l- 



and the numerical factors ai^ are: 



(4.7) 



03 



327r\^ V2 



ai = - 



8V27rH2 



For two-dimensional systems the result is, with loga- 
rithmic accuracy: 



T 
hD 



(l + 3A2)lln(^ 



+ (1 + 3A2) ^C(3)Tr 



(4.8) 



where C(3) ~ 1.202. The first term on the right hand 
side extends to higher temperatures the logarithmic be- 
haviour known in the diffusive limit (the upper cutoff is 
of the order of the Fermi energy Ep and not ?i/t); the 
second term becomes relevant in the quasi ballistic limit. 
In the diffusive limit our results are the same as those 
obtained in Ref. 12 by explicit thermodynamic calcula- 
tion. 



V. DERIVATION OF THE KINETIC EQUATION 

This section is devoted to the derivation of the local ki- 
netic equation. We first introduce the Eilenberger equa- 
tion and some basic notation. Then we perform a (gen- 
eralized) gauge transformation: this is the crucial step 
that enable us to obtain the local description. Then we 
introduce the bosonic degrees of freedom and derive the 
collision integrals. 



A. Eilenberger equation 

Our starting point for the derivation of the kinetic 
equation is the same as in Ref. 17 and we briefly summa- 
rize it here. 

The interaction with small momentum and energy 
transfer in the singlet channel [the triplet channel will 
be discussed in section VH] is decoupled using the two 
Hubbard-Stratonovich fields 4>±(t, r, n). For the purpose 
of the one loop approximation we will employ, these fields 
can be considered as Gaussian with the propagators 



((<^+(l).^_(2))) = --D«(l,2), 

((0_(1)0_(2)))=O. 



(5.1) 



where ((...)) means averaging over the fields (j)±. We used 
the short hand notation 



drii 



(i*) = {ruTii), j di* = j dri j ■ 
where i = 1,2, . . . and is the total solid angle. 



drii 



(5.2) 
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We introduce the disorder averaged Green's function of 
the electron in the field (j)± in its matrix form in Keldysh 
space: 

such that its average over the fluctuating field ^± gives 
the usual expressions for the physical propagators: 

((G«(l, 2))) = -i9{h - i2)(V(l)V'^(2) + V^(2)V(1)), 

((G^(l,2))) = z0{h - <i)(V'(l)V'^(2) + ^-^(2)^(1)), 
((G^(l, 2))) = -^(V(1)V^(2) - i^Hm^)}, 



((G^(l,2)))=0. 



(5.4) 



Here 9{t) is the Heaviside step function, ip^ , tjj are the 
fermionic creation/annihilation operators in the Heisen- 
berg representation and quantum mechanical averaging 
(. . .) is performed with an arbitrary distribution function 
to be found from the solution of the kinetic equation. 

For the disorder averaged Green's function, the semi- 
classical approximation is obtained by integrating the 
Wigner transform of G{1, 2\^) over the distance from the 
Fermi surface: 



G{tut2,p,R) = d'r e'^'-G(l,2 



(5.5) 



ri - r^; 



R 



-Xri + r2); 



P = p--[A {ti,R) + A{t2,R)]; 



g{ti,t2,n. 



OO 



d£, G ti,t2,n 



Pf 



Vf 



(5.6) 

where A is the vector potential of an external electromag- 
netic field, Pf the Fermi momentum and Vf the Fermi 

velocity. The dynamics of the semiclassical Green's func- 
tion g in the matrix form is governed by the Eilenberger 
equation^^: 



dt + v - V + u>c 



dn 



g ° Strg 



(5.7) 

where v = vfu, the action of the "bare" collision integral 
on any function a{n) is defined as 



Sti-a 



, , f drii „ , N / \ 
J "Q^^*'^v'^>^i)a("i) > 



Str(ni,n2) 



r(^i2) 



S (nin2) 



dn2 

t{9i2) 



(5.6 



and 9i2 = nin2 [for the short-range impurity t{9) does 
not depend on 9; however, the formulas derived here will 



be valid for the arbitrary impurity scattering] . The time 
convolution of two matrices a(ti,<2) and h{ti,t2) is: 



aob = Jdts d{ti,t3)b{t3,t2); 



d ° b = dob — bod. 

(5.9) 

Defining the commutator between a matrix c{t, r, n) and 
9 as 

[c; g] = c{ti,r, n)g{ti,t2,n, r) - g{ti,t2, n, r)c{t2,r, n) 

(5.10) 

the covariant derivatives in Eq. (5.7) are given by: 



dt'g = dtj + dtj + i[(p;g] 



Vg = Vg + i 



A; 5 



(5.11a) 
(5.11b) 



with A = Al/f and — (fii-K- Here Ik denotes the unit 
matrix in Keldysh space and is the scalar potential for 
an external electromagnetic field such that: 



eE = -V(p + dtA 



°B = -cV X A 



The vector = eB/{mc) has the magnitude of the cy- 
clotron frequency and the direction of the magnetic field 
B. Finally, ^ is the matrix in the Keldysh space: 



0- (t>+ 



(5.12) 



K 



The matrix Green's function g is subject to the follow- 
ing constraints: 



,g(n, r)og{n, r) = S{ti - t2)iK 
dt Tr g{t,t,n,r) = 0. 



(5.13a) 
(5.13b) 



In thermal equilibrium, the relation 

g'^itiM) = [g^^o n-no g^] {tut2) 



n{ti,t2) 



de 
2^' 



,ie{t2-ti) 



n(e); 



(5.14) 



n(e) = 1 - 2/F(e) = 2tanh 



2T 



must hold independently of the form of the spectral func- 
tions g^'^. 

In what follows we will assume that there is no mag- 
netic field, B = and u)c = 0, but no gauge choice is 
made: although one could set A = by a gauge trans- 
formation, both the scalar and vector external potentials 
are left arbitrary in order to keep track of the gauge in- 
variance of the equations. 

As for the propagators defined by Eq. (5.1), they sat- 
isfy matrix Dyson's equation: 



D{1, 2) = £)o(l, 2) + jdS jdA Do{l, 3)n(3, 4)£)(4, 2) 
FPi9i2)Siri2) 



Do(l,2) 



V{ri2) + 



S{ti2)tK; 

(5.15) 
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where V{r) is the long range part of the interaction (for 
the Coulomb interaction V{r) = e^/|r|), 6*12 = «a«^2, 
ri2 = Ti— r2, ti2 = ti — t2- The matrix propagator is 
denoted by D and n is the matrix polarization operator. 
They have a structure similar to the Green's function 
one: 



^ - ' 



n 



K 







(5.16) 



K 



The polarization operators are given by variational 
derivatives of the solutions to the Eilenberger equation 
(5.7): 



n^(i,2) = n^(2,i) = zy 



512 



2 S(j)+{t2,r2,n2) 

(5.17a) 



n^(l,2) = 



TTi/ 



(5.17b) 



(5.18) 



S(l)-{t2,r2,n2) 
where 

S12 = ^dS (nina) S{ri - r2)5{ti - t2), 
with being the total solid angle. 

B. The gauge transformation 



With the Eilenberger equation (5.7) at hand, one could 
proceed as in Rcf. 17 in order to derive an equation for 
the distribution function. However the resulting inelas- 
tic part of the collision integral, expressed in terms of 
the electron distribution function only, is non-local and 
the evaluation of e.g. the thermal conductivity would 
require the time and spatial gradient expansion of this 
term in the spirit of Eq. (2.15). As we already discussed, 
such a route makes the energy conservation in the kinetic 
equation obscure. Here we follow a different approach, 
inspired by the following considerations-^^: if the fluctu- 
ating fields were uniform, they could be eliminated from 
Eq. (5.7) by a gauge transformation 



(5.19) 



In other words, the position-independent fluctuations of 
the <^ fields define the time dependent position of the en- 
ergy levels but the occupation numbers for such levels 
do not change. Therefore such fluctuations affect nei- 
ther the electric transport nor the electron contribution 
to the thermal transport in the system. Moreover, if the 
path of the electron were a straight line, all the smooth 
fluctuating flclds could still be eliminated in the cikonal 
approximation and once again they should not affect the 
electronic contribution to the transport. To get rid of 
such spurious contributions we will employ the gauge 
transformation described below. 



We introduce a new matrix field K{t, n, r): 



K = 



(5.20) 



which is a functional of the field cj) and is used to perform 
the "generalized" gauge transformation 



g-iK(ti ,n,r) - giK(t2 ,n,r) 



(5.21) 



This transformation is unitary and as such preserves the 
constraints (5.13). As we will sec, it will lead us to the 
local kinetic equations. Applying the transformation to 
the Eilenberger equation (5.7), we obtain 



dt + v-V 

where 

■^4, ■ 
St^.g 



{dt + v-V)K-cf>,g 



drii 



9 ° St^.gJ , 
(5.22) 



{ti,t2,n) = J -^St^(n,ni) 



The "bare" impurity collision integral and the derivatives 
are defined respectively in Eq. (5.8) and (5.11). We sup- 
pressed the r argument which is the same in all functions. 

We will look for a perturbative solution to the Eilen- 
berger equation in the form (5.22) in the first loop ap- 
proximation; to do so it will suffice to retain only the 
terms at most quadratic in the K fields in the collision 
integral. At the lowest order, g has the form: 



9 = 



_ fS{h - t2) 



K 
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-S{h - t2) 



(5.23) 



We will require that such form is preserved even in the 
first and second order in i.e. the corrections to the 
spectrum (described by g^'^) are indeed eliminated by 
the gauge transformation. 

At the linear order, the retarded, advanced and 'Z' 
components of Eq. (5.22) vanish if K- satisfies the equa- 
tion 



(5.24) 



The solution to the intcgro-differential equation (5.24) 
can be written in terms of the diffuson propagator 
i^^{ti-,t2',ni,n2;ri,r2) - the retarded solution to the 
classical kinetic equation 



(dt, + VI ■ Vr, - St,) Ca = S12, 



(5.25) 



where 612 is defined in Eq. (5.18). Using Eq. (5.25), we 
find 



K_{1) 



d2£5(l,2)</._(2), £3(1,2) =£5(2,1), 

(5.26) 
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where we use the shorthand notation (5.2). In an opera- 
tor notation, Eq. (5.26) can be rewritten as 



(5.26') 



To simphfy further manipulations, we introduce the 
following function of three angular variables 
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(5.27) 

This function is related to the impurity collision integral 
(5.8) by 



dn2 
" dni 



ni;n2 

"-3 

ni;n2 

"•3 

ni;n2 



[St^] (ni,n3); 
- [Sir] (n2,n3); 
0. 



(5.28) 



Denoting with 6g^ the first order correction to , the 
Keldysh component of the Eilenberger equations at the 
linear order is: 



i(dt + v- V - St, 



K 



(5.29) 



dt+v-V -StrjK+-(t)+\ 



QK. 



with [cf. Eq. (5.11)] 

dtSg'' = (dt, + dt,)Sg'' + Sg""] 
VSg"" = VSg'' + i[A;Sg'']. 

Here Q is a local in space operator 

QK- = jdta j ^^Q{ti,t2,n;t3,nr,r)K_{t3,ni,r) 
with the kernel 



Qiti,t2,ni;t3,n2) = - J ^7 



n2 



yni;n3j (5.30) 
X [Pg^ {ti,t3,ni)g^ {t3,t2,n3) + (ni ^ ns)] . 

We suppressed the spectator argument r, which is the 
same in each term of the equation, and the last term 
means that one must add terms like the ones shown 
but with the angular arguments of the Green's function 
switched. The principal value sign P in Eq. (5.30) means 
that the divergent at ti — > t2 part of the product of the 
Green's functions 



g^{ti,t2,n,r] 



2i 



ti^t2 7r(ti - 12) 



regular (5.31) 



must be excluded: 
Pg''{ti,t3)g''{t3,t2) 

= 9''{tut3)g''{t3,t2) - 45(tl - t3)6{t3 - t2), 

or, equivalently, 

P9''itut3)g''{t3,t2) 

= ^ V g'^{ti,t3 + aiO)g''{t3 + <7iO,t2). 

It is worth noticing that all non-equilibrium effects con- 
tribute to the regular part in Eq. (5.31) but not to the 
singular part - the states deep into the Fermi sea, which 
are not perturbed, contribute to it. 

In order to solve Eq. (5.29) we define a new field K- 
by the relation 



k_{t,n,r) = (idt) MK_ , 



(5.32) 



where the operator M will be shown to be related to cer- 
tain products of the Green's functions g^ , see Eq. (5.38). 
The operator M is Hermitian and local in space but not 
in the momentum direction and time. We used again an 
operator notation 



MK 



j dti j ^^M{t,n;ti,ni;r)K_{t-i,ni,r). 

(5.33) 

We require K+ to satisfy: 

(5.34) 



(dt+v-V- Str^ K+ = (P+- 2k. 

whose solution is: 

K+= t^(j)+-2t^k_. 



(5.35) 



The operator notation here is the same as in Eq. (5.26). 

The next task is to choose the "best" form for the op- 
erator M to maximally simplify any further perturbative 
expansion. Writing 5g^ = 6g^ + 6g^ , we obtain the 
following equations: 



iLSg^ 
iL5g^ 



QK_+2 

dn2dn3 



K-,g 



K 



I 



^2 



-7 



"3 

n;n2 



(5.36a) 

(5.36b) 



X [i^+(n3),(ff^(n2)-5^(n))], 



where L= (dt + v-V - St^) . 

We note that the right hand side of Eq. (5.36b) van- 
ishes in equilibrium, since g^ = {g^)n- Therefore 5g^ 
also vanishes in equilibrium and cannot contribute to 

equilibrium properties such as the specific heat. More- 
over, even in non-equilibrium 6g^ {ti,ti,n,r) = 0, be- 
cause the right hand side of Eq. (5.36b) vanishes, see the 
remark after Eq. (5.31). It means that 5g^ does not 
contribute to the electron density or current. 
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We are now ready to choose the operator M. We re- 
quire that Sg^ {ti,ti,n,r) = 0, i.e. Sg^ also does not 
contribute to the electron density or current. It means 
the the right-hand side of Eq. (5.36a) must vanish for 
ti = t2 for any field Imposing this requirement, we 
obtain 



M{ti,t2,n,fi,r) = -Q{ti,ti,n;t2,n;r) 



Together with Eq. (5.30), it yields 



n2 



(5.37) 



(5.38) 



X [Pg^ {ti,t2,ni)g^ {t2,ti,n3) + (m ^ ns)] . 

Expression (5.27) for the vertex 7 enables us to establish 
the following properties of kernel Eq. (5.38) 

M{ti,ni;t2,n2) = M{ti,n2\t2,ni) = M{t2,n2;ti,ni), 

(5.38') 

i.e. the operator M is Hermitian. 

It is instructive to find M in the thermal equilibrium. 
Using Eq. (5.14) and the fact that due to the choice (5.24) 

the retarded and advanced components of g are still given 
by Eq. (5.23), we find from Eqs. (5.28) and (5.38): 

_^e^'^(*^-*i)Meg(a;;ni,n2) 



and we use the notation (5.2) throughout this subsection. 
The result for the Keldysh component is 



n^(l,2) = 2w C^MjC^ (1,2) 



(5.43b) 



The actions of the operators M and D are defined in 
Eqs. (5.26) and (5.33). 

It is easy to check that the fliictuation-dissipation re- 
lation between the polarization operators holds in the 
thermal equilibrium. As follows from Eqs. (5.25) and 
(5.26) 



(5.44) 



We perform the time Fourier transform for all the propa- 
gators and the polarization operators in thermodynamic 
equilibrium 



Ail 



,iuj(t2-ii) 



A(w;r,2*), 



(5.45) 



where we use the short hand notation (5.2). Substituting 
Eqs. (5.39) and (5.44) into Eq. (5.43b), we obtain that 
in equilibrium; 



Meg(a;;ni,n2) = — wcoth 



f- 

V2r 



Str 



(ni,n2). (5.39) nf,(a;; 1*, 2*) = n«(a;; 1* 2*) - n^(a;; 1* 2*) 



Equation (5.39) will be useful for checking the fluctuation 
dissipation theorem. 



UJ 

coth — -. 
2T 

(5.46) 

and with the help of Eq. (5.15), we derive the fluctuation- 
dissipation relation 



C. Polarization operators and propagators 

The knowledge of the linear order corrections to the 

Green's function permits the calculation of the polariza- 
tion operators as variational derivatives of the original 
Green's functions (i.e. before the gauge transformation) 
in the limit t2 ii, see Eq. (5.17). At linear order the 
corrections to the original Green's functions are given by 
the relations [cf. Eq. (5.21)]: 

dg"" ^ Sg"" -i[K+,g'^]-2K_5{H-t2), (5.40a) 

5g^ 2K_6{ti - t2). (5.40b) 

By construction of the previous subsection: 



lim 5g^ {ti,t2,n,r) = Q 

(2^*1 



(5.41) 



and using Eq. (5.31): 



lim —i 



[K+,9 



--dtK+{t,n,r). 

TT 



Substituting these results into Eq. (5.17) and using 
Eqs. (5.26), (5.32) and (5.35), we obtain: 



U^{2,1) = u [6,2 -dt,C<^ {1,2)], 



(5.43a) 



D^^ico; 1*,2*) = [O'^iw; 1*,2*) - D^{u; 1*2*)] coth^. 

(5.47) 

Given the expressions for the polarization operators 
obtained above, we can solve the Dyson equation (5.15) 
and obtain the explicit expressions for the interaction 
propagators. In the operator notation: 



= — 



1 



l + F- dtFCs 
1 



z/D^ = -F— 

1 + F + dtCsp 

= 2iub^jo.3M£3f)A 



(5.48a) 
(5.48b) 
(5.48c) 



where the action of the operator F on any function 
a{t, n, r) is defined by 



(5.42) [Fa]{t;n,r) = J ^ 



FP{nrii)a{t, ni,r) 



+ 



J driuVii — ri)a{t,ni,ri) 



(5.49) 



see also the text after Eq. (5.15). 

To find the propagators for the fields K± given in 
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Eqs. (5.26') and (5.35), defined as: 

((if+(l)i^+(2))) = ^/C^(l,2), 

((i^+(l)i^_(2))) = ^/C«(l,2), 

{{K_{l)K+{2)))='-lC^{l,2), 
{{K_{l)K_{2)))=Q 



(5.50) 



we use Eqs. (5.1) and (5.48) and obtain for the retarded 
and advanced propagators 

= ab^a- = ab^a, (5.51a) 

whereas the result for the Keldysh propagator is 

t^=-t9b'^& 



+ 2i 



CB{dt)-^Mab^CB - jCBb^&M{dt)-'^C<' 

(5.51b) 



The fluctuation-dissipation relation between the D pro- 
pagators in Eq. (5.47), the equilibrium form for M in 
Eq. (5.39) and the identity (5.44) enable us to verify the 
fluctuation-dissipation relation for the /C propagators: 



/C^(a;; 1* 2*) = [/C^(a;; 1* 2*) - ]C^{co; 1*, 2*)] coth 



2T 
(5.52) 



D. Additional bosonic fields 



Equation (5.51b) together with Eqs. (5.51a) and (5.38) 

permits to express the Keldysh propagator JC^ in terms 
of the electron distribution function. This relation, how- 
ever, would be nonlocal on the spatial scale much larger 
than the temperature length 



i7 



hv p I Tit 



(5.53) 



recall the discussion of Sec. II B. Indeed the collision 
integral and all the physical quantities will be given by 
integrals of the type 



j duf{uj)}C°'{io), a = R,A, K, 



where the function /(w) depends on its argument on the 
characteristic scale of T. A retarded function is an an- 
alytic function of u at Imu > 0, which means that for 
a = R the integral will be determined only by the sin- 
gularities of /(lu), i.e. ~ IC^{lu = iT). This im- 
mediately restricts the spatial scales to Lt- The same 
argument applies to the advanced case, because of ana- 
lyticity at Imw < 0. The function IC^ (ui), however, is not 
analytic. Moreover, according to Eq. (5.51b) it contains 



overlapping singularities of the retarded and advanced 
propagators. It means that the characteristic frequencies 
entering T^"^ are determined by the poles of the propaga- 
tor rather than by the width of the function /, i.e. the 
spatial scale may by far exceed Lt and any expression of 
the type is thus non-local. 

To overcome this difficulty, the standard parametriza- 
tion of the Keldysh function = D^oN - NoD^ is 
usually introduced and the kinetic equation for the dis- 
tribution function N is then derived. All the non-locality 
in the problem is then contained in the partial solution of 
the kinetic equation, compare with Eq. (2.14), whereas 
the kinetic equation itself is local. 

In what follows, we adopt this program in a slightly 
modified form. We introduce a new retarded propagator 
£''(1,2), compare with Eq. (5.25) 



iH,_iSdi,,-iVi) -Str £" = (512, 

CO (5-54) 



He-h{iO,q) = vq- 



1 + F 



and its advanced counterpart £''(1,2) = £''(2,1). The 
multiplications in Eq. (5.54) are to be understood in the 
operator sense and the action of the operator F on a 
function a{t,n,r) is defined by Eq. (5.49). 

For F = 0, C (£'') reduces to the usual diffuson £f 
(£^). Physically, C describes the spectrum of the prop- 
agating electron-hole pair and the operator in the left 
hand side of Eq. (5.54) corresponds to the kinetic equa- 
tion for the collective mode in the Fermi liquid theory^^. 
The operator Hc_\i{uj, q) can be interpreted as a "Hamil- 
tonian" [see also App. C] of the interacting electron-hole 
pair. 

In terms of £'' and £^, Eqs. (5.51) acquire the form 



uK,^ = - -CP (dt) 



K 



(5.55a) 



{dt)-^. (5.55b) 



Let us introduce two bosonic "distribution functions" 
(the density matrices to be more precise), and N^, 
that satisfy the following equations: 



{c^y^N^^ + ^f3{C3)- 



= 2M, 
= 2M. 



(5.56a) 
(5.56b) 



The operator M is defined in Eqs. (5.33) and (5.38) and, 
in a more explicit notation, the action of the operators 
NP'9 on any function a{t, n, r) is to be understood as 



(1) = jd2AfP'S (1,2) a{2), 



where the shorthand notation (5.2) is used. Note that 
from the above equations it follows that the bosonic func- 
tions AfP'9 are symmetric 



AfP^^l, 2) =AfP'^ {2,1). 



(5.57) 
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Equations (5.56) enable us to rewrite Eq. (5.55b) as form the time Wigner transform: 



uJC" =-i {dtY 

+ i{dty 



+J(;f9jr9 
CPMP +MPCP 



(5.58) 



This expression is local in the above discussed sense and 
will be used in the construction of the conserved energy 
current. Obtaining the local expression, however, re- 
quired the introduction of two additional bosonic distri- 
bution functions: M^, describing the interacting electron- 
hole pairs and the ghost field distribution Af^, subtract- 
ing the contribution of the electron-hole pairs in the ab- 
sence of interactions. 

Closing this subsection, let us rewrite Eqs. (5.56) in a 
form resembling the kinetic equation of Sec. II. We sub- 
stitute Eqs. (5.25) and (5.54) in Eqs. (5.56) and obtain 



dt + v-V; Af3] = St {A/'^.g^} (5.59a) 



iH,.^{idt,,-iVi); MP =St {A^^5^},(5.59b) 



where the collision integrals are 



St 



{AA",c/^} = 2{st^;Ar"} +2M; (5.59c) 



for a = g,p. They depend on via M and we used the 
notation 

^^^A;B^ = ^{AB + BA); A;B = AB - BA. (5.60) 



We perform the time and space Wigner transforms of 
Eqs. (5.59) to introduce the bosonic distribution func- 
tions NS'P 



" giq-(i'i-i-2) 

{2rrf 



X U) 



[2NS'P {u, q; m, n2; t, r) + QdS {nrn2) ] , 



(5.61) 



where t = r ^ 

(5.57) translates into the condition 



^^^-ira. The symmetry relation 



N^'P (a;,q;ni,n2) 

= - [NS'P (-w, -q; 712, «i) + (nlna) 



(5.62) 



The physical meaning of this relation is boson statistics: 
at w > 0, N^'P corresponds to the occupation numbers 
entering the probability of the absorption of the bosons 
whereas the iv < part describes the boson emission. 

The fermionic distribution function / is obtained in 
two steps: (a) we introduce the gauge invariant Green's 
function g (see also the next subsection) and (b) we per- 



g {ti,t2,n,r) = exp 
g{ti,t2,n,r) = 2 



fde 



-t j dtip{t,r)j g{ti,t2,n,r) 
(5.63a) 

(5.63b) 



Performing such Wigner transforms on Eqs. (5.56) and 

(5.38), we find 



St {N^f} 

(5.64a) 



1 



l + F 
+i 



■,d,NP 



[s{u;,q);VNP} 



H,.^{uj,q),NP 



= St {NPJ}, (5.64b) 



where the collective mode velocity operator is 

dHMu;,q) d f F \ 

s{w, q) = = V + LJ 



dq 



dq V 1 + F 



(5.65) 



In the left hand side of equation (5.64b) we limited our- 
selves to the leading Poisson brackets [the equation be- 
comes exact for the short range interaction, since dqF = 
0, and in the unitary limit F oo]. However no Poisson 
brackets arise in the right hand sides of Eqs. (5.64) - this 
is a consequence of the locality of the kinetic equations. 

The right hand sides of Eqs. (5.64) describe the de- 
cay of the electron-hole pair into the electron and the 
hole moving into different directions. To write down the 
expression for this collision term, it is convenient to in- 
troduce the following object: 

'^iifw 9! = ('^' 9; ni,nj;t, r) 

x{/(e,fife;i,»') - /(e-t^jTifc ;*,»•)} (5.66) 

+Q.d5{n^j)^f{e,ni;t,r) [1 - /(e - a;,nfc;i,r)] |. 

It is easy to see that T^''' = in the thermal equilibrium 
Eq. (3.1). 

In terms of this object and the vertex (5.27) we have 
Si' {A^«■^ /} (a;, q;nun2;t,r) = -jdej 



(5.67) 



where we suppressed the spectator arguments t, r and 
q in the right-hand side of the equation. Deriving 
Eq. (5.67) we used Eqs. (5.28), (5.38') and the property 

Jde[f{e)-f{e-co)]=-uj. 
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To understand the physical meaning of the processes 
described by the coUision integral (5.67), we use the ex- 
plicit form of the vertex 7 [Eq. (5.27)] for the isotropic 
impurity scattering t{9i2) = t. Then, the collision inte- 
grals can be decomposed into the sum of two contribu- 
tions 

StiN^'", /} = S^iN^'", /} + StliN^'", /} 

The first term in the right-hand side can be obtained from 
a simple coimting of the probabilities of the processes 
depicted on Fig. 2(a-b) 

StliiN^'P,/} (c.;ni,n2) 

= 7 / ^ { "^12-32 (e, w) Tf ^'^la (e, w, q) } . 

The second term in the right-hand side originates from 
the interference of two scattering processes, see Fig. 2(c). 
That is why it produces contributions to N which are not 
diagonal in the momentum directions: 

E. The collision integral for electrons 

With the bosonic propagators /C at hand, we can pro- 
ceed with the calculation and include the second order 
in the fluctuating fields K± contributions to the colli- 
sion term of the Eilenberger equation (5.22). With the 
fluctuating fields K± given by Eqs. (5.26') and (5.35), 
Eilenberger equation becomes: 



1 



St^g ■ 



(strK+ - 2^_) Ik + iStrK^a'^K 



(5.68) 



where we use the notation (5.11) for the derivatives, St^ 
was defined after Eq. (5.22) and K- in Eq. (5.32), and 
is the Pauli matrix. 

We expand the right hand side of Eq. (5.68) up to 
the second order in K, see Fig. 3a; then we average it 
to obtain Fig. 3b. The resulting second order contribu- 
tions can have two different origins: (1) they may arise 
from the expansion of the exponentials truncated at the 
second order, term Sti on Fig. 3c, or (2) they are ob- 
tained as products between the linear correction 6g^ of 
Eqs. (5.36a) and the first order expansion of the expo- 
nentials, term St2 on Fig. 3c. 

The Eilenberger equation for the averaged Green func- 
tion takes the form 



dt + v-V 



[dt + v- v) 



a) 



b) 



K K 



(dt + v- v) 



C) 



= 

St St^ Sti 



St2 



- w ■ V - St, 



FIG. 3: Schematic representation of the averaging over the 
fluctuation fields K. a) Expansion of the Eilenberger equation 

(5.68) before the averaging; b) The equation for the Green 
function averaged over K, see Eq. (5.69); c) The contribtu- 
tions to the total collision integral in the first loop approxi- 
mation. 



where St contains both zero and second order contribu- 
tions. We find [see Appendix D for details on the can- 
cellation of the second order corrections in the R,A,Z 
sectors] : 



St K, 50^,5^} 



St {g^,NP,N3} 




(5.70) 



where g^ = —g^ = 5{ti — 12)- It means that the matrix 
Green's function of the form (5.23) is still the solution 
of the Eilenberger equation - the main gain of the gauge 
transformation (5.21) - provided that the kinetic equa- 
tion for the Keldysh component is satisfied: accordingly, 
we concentrate on this component only. 

Performing the gauge transformation (5.63a) on the 
Keldysh component of the Eilenberger equation, we ar- 
rive at the explicitly gauge-invariant form of the kinetic 
equation: 



g = St{g,NP,N^ 



g = Si{g,NP,Na} 



(5.69) 



dt.+dt^+v-V + iv-f dieE{i,r) 

(5.71) 

In the collision integral St{g,NP,N^}, we indicate 
with Sti{g, A'''', N^} the contributions of the type (1) and 
with St2{g,NP,N3} those of type (2) [this separation 
has no particular physical meaning - it is just a matter 
of practicality in the calculations; we will return on the 
physical aspects when analyzing the conservation laws in 
the next subsection]: 

St{g, N", N^} = St,5-4Sti{5, A^^, NS}-^St2{g, N", N^}. 

(5.72) 

where St^ is defined in Eq. (5.8). The numerical factors 
in front of the last two terms are introduced to facilitate 
the transformation to the canonical form of the kinetic 
equation in further sections. 
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The expression for Sti written in terms of the K prop- 
agators (5.50) and the 7-vertex (5.27) is 



H (^i,^.;ni) = -^/^7l\ 



fields K± with the propagators (5.50) and 

St-{ti,t2;n,r) = ^ f dt^ f ^^K-{t3,ni,r) 



(5.73) 



X {g{h, t2, na) [jC§,{h,t2) - iCg{h, 12) 

+ Jdt3g{ti,t3,ni)g{t3,t2,n2) [/C^2(*3,t2) - /C^i(i3, t2)] 

-9{tl,t3,n2)g{t3,t2,ni) [/Cij(ti, ig) - /C^Kii, fg)] }■ 



6Q{ti,t2,n;t3,ni;r) 



(5.76b) 



+ {ti,t2;n,r) / dt^Q{t^,h,n]t3,ni;r) 



We introduced the short hand notations 



St+(ti,t2«') = T / — — 7 



ii'+(n3); (<^ff(n2) - (55((n)) (ti,t2). 



(5.76c) 



}Cij{ti,t2) = IC{ti,ni;t2,nj), 



(5.74) 



and 



The commutator was defined in Eq. (5.10) and the kernel 

5Q is the first variation of the operator (5.30) with re- 
spect to the Keldysh component of the electron Green's 
function 



IC^{t„t2)^2IC^{t„t2)-IC^{h,h)-K:^{t2,t2). SQ{t„t2,n;t3,n,;r) = \J ^7 (^^^ J 



We omitted the variable r, which always appears in the 
distribution function as g{ti,t2,n,r) and in the propa- 
gators as IC{ti,ni,r;t2,n2,r). The dependence on the 
electron distribution function g is explicit, whereas the 
dependence on the bosonic distribution functions is hid- 
den into the propagators - see Eq. (5.58). For reasons 
which will become clear in the next subsection, we split 
Eq. (5.73) into two parts 



{ 



Sg+{ti,t3, n, r)g{t3,t2, n2, r) (5.76d) 
+g{ti,t3,n,r)6g+{t3,t2,n2,r) + {n2 ^ n) 



dn2dn3 3 



(5.75a) 
(5.75b) 



Sti = Sti + Sti ; 
[St'i"] (ti,f2;ni) = -^ J 

X { 25 (t 1 , i2 , n2 ) [^^ (t 1 , i2 ) - ^fi (^1,^2) 

+Jdt3 [g{ti,t3,ni)g{t3,t2,n2)+ g{ti,t3,n2)git3,t2,ni)] 

X [/C4(i3,i2)-/Ci(ii,t3)-/C^i(t3,i2)+/Cf3(ti,i3)] }. 



Sti 



, , i f , f dn2dn3 , . . 

(ii,t2;ni) = -- Jdt3 J — ^7?2 (5.75c) 

x{ [9ih,t3,ni)g{t3,t2,n2) - g{ti,t3,n2)g{t3,t2,ni)] 

X [lCi2{t3,t2)+ICi3ih,t3)-ICi^{t3,t2)-ICUh,t3)] }. 



As for St2, it is convenient to separate it in two parts, 
depending on which field, or we retain in the 
expansion: 



St2 = ((St+)) + ((St_ 



(5.76a) 



where ((. . .)) stands for the averaging over the fluctuating 



Finally, the functions Sg± and 5g = Sg^ + 5g_ are 
obtained by solving Eqs. (5.36a) and (5.3Gb) [after the 
transformation (5.63a)]; with the help of Eq. (5.37): 

Sg+{tut2;n,r) = 'i J J ^^73^2 (5-77) 
£^ (ti,n,r;3) ^^+(^3, 714, ^3) - K^{t3 - ti2,n4,r3) 
X g^{t3,t3 - ti2; n2, rs) - g^{t3,t3 - ti2; n3, r3) 



5g-{h,t2;n,r) = -i J d3 / '^^4 ^ ^ (5.78) 
£^{ti,n,r;3) K_{t4,n4,r3) Q(<3,<3 - ii2, ".3; ^4, ^4; r3). 

-yfi'^(*3,t3 -ti2;n3,r3) j dt^Q{t5,t5,n3;ti,ni;r3) 

t3—tl2 

where ti2 =ti — t2 and the notation (5.2) is used^^. 

For future use, we note following properties [see also 
the discussion following Eq. (5.36)]: 

5g-{t,t) = 5g+{t,t) =^'': j dn5Q{ti,t2,n;t3,ni;r) = 0, 

(5.79) 

as one can see from Eqs. (5.77)-(5.78) and from the defi- 
nitions (5.76d) and (5.27). 

The canonical form of the kinetic equation is obtained 
by performing the time Wigner transform (5.63b) of both 
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sides of Eq. (5.71). It is clear from the structure of the 
collision integrals that this procedure will lead to the ap- 
pearance of Poisson brackets in the right hand side of the 
kinetic equation. We choose, however, another route: we 
will prove the existence of the conservation laws before 
the Wigner transform. It will enable us to argue that 
these Poisson brackets (in our formulation of the kinetic 
equation) give only small contributions, which can be ne- 
glected within the accuracy of the kinetic equation. 



F. Conservation laws 

The derivation of the conservation laws is based on 
the following properties of the collision integrals of the 
previous subsection 



(St-r5)r 
(St_)„ 

(Stf)„ 



= 0; 
0; 
0. 



(5.80a) 
(5.80b) 
(5.80c) 



The physical meaning of conditions (5.80) is that the 
corresponding terms in the collision integral conserve the 
number of particles within the energy shell. Equations 
(5.80) follow immediately from definitions (5.8), (5.75b), 
(5.76b), (5.76d) and (5.79). 

The two remaining contributions to the collision inte- 
gral have the properties 



lim St+(ii,t2) =0; 



lim Sti (ii,t2) = 0, 



(5.81a) 
(5.81b) 



and 



lim (dt, - dt,) St+{tut2) = 0; 

tl— »t2 



lim {dt,-dt,) Sti itut2) 

tl— ►12 



(5.82a) 
(5.82b) 

jdts [dttC^] ^3 {ti,t3)[g{ti,t3,ni)g{t3,ti,n2) 

+g{ti,t3,n2)gih,ti,ni)] - (ni n2)|, 

where wc use the notation (5.74) and the vertex is de- 
fined in Eq. (5.28). Equations (5.81a) and (5.82a) imme- 
diately follow from the definition (5.76c) and the condi- 
tion (5.79). Derivations of Eqs. (5.81b) and (5.82b) are 
presented in Appendix E. 

Expressions (5.81) mean that while not conserving the 
number of particles for a given energy shell, the terms 

^ ^ in 

St_|_ and Sti conserve the total number of particles for a 
given direction (small angle inelastic scattering). Equa- 
tion (5.82a) means that the inelastic St+ term conserves 
not only the number of particles but also the energy for a 



given direction. Equation (5.82b) means Sti term does 
not conserve the energy for a given direction, thus de- 
scribing the energy exchange between the quasiparticles 
and the electron-hole pairs discussed in Sec. VC. 

The possibility to find the conserved energy current is 
based on a certain relation between Eq. (5.82b) and the 
collision integral for the electron-hole pairs - we now turn 
to the discussion of this relation. 

We substitute Eqs. (5.55a) and (5.58) into Eq. (5.82b) 
and average the result over n; then using the analytical 
property 



,.) = / 



2^ 



IC^{u) = 



(5.83) 



together with £''(1,2) = £''(2,1), £^(1,2) = £^(2,1), 
J\fP'3{l, 2) = N'P'3{2, 1) and Eq. (5.38), we find: 



hm {dt,-dt,){Sti % 

tl— >t2 



(5.84) 



= TV„ 

+— TV„ 



a (X/'sst^ + m) 

where St^ is defined in Eq. (5.8) and Tr„ acts as 



Tr„A 



(i,r) = j ^A{t,n,r;t,n,r). (5.85) 



The corresponding traces of the collision integrals for 
the electron-hole pairs - Eq. (5.59c)'*° - are: 



Tr^ 



Tr^ 



St" {Af<^,g} 



Tr, 



+2C^ Af^Str + M 



&St' {AfP,g} 



Tr„ 



CP, St, 



(5.86a) 



AfP 



+2CP { AfPStr + M 



(5.86b) 



Comparing Eqs. (5.84) and (5.86) and using once again 
Eqs. (5.25) and (5.54), we arrive to the desired relation 
between the collision integrals: 



2i7rz/ lim {dt, - dt,) {St {g,J\fP ,Af<^})r 



(5.87) 



+TV„ 
= Tr, 



£PSt'{AfP,g} 



-Tr. 



£'^St"{Af'^,g} 



H,.^{idu-iV),tP 



MP 



The left hand side of Eq. (5.87) is the quantum coim- 
terpart of the relation (2.3b) and Eqs. (5.80)-(5.81) are 
related to Eq. (2.3a); we will derive the expressions for 
the electric and energy currents in the spirit of our dis- 
cussion in Sec. II A. 
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Wc begin with the conservation of electric charge. Ac- 
cording to Eqs. (5.4) and (5.6) the charge density is given 

by 

P{t,r) = — hp {g{ti,t2,n,r))r,. (5.88) 

Taking the hmit ti ^ t2 ^ t in both sides of Eq. (5.71) 
and using Eqs. (5.72), (5.80) and (5.81), we obtain the 
continuity equation 

5tp+V-j = 0, (5.89) 

where 

j{t,r) = — hm {ng{ti,t2,n,r))n, (5.90) 

Z tl — *r2 — *t 

compare with Eqs. (2.5) ~ (2.7). 

Having found the usual equation for the electric cur- 
rents, we turn to the energy conservation. Acting with 
the operator {dt^ — dt^) on both sides of Eq. (5.71) and 
introducing the quantities 



where the energy densities Up^g and currents ^ are de- 
fined as 

Ug{t,r) = iTr„£s J^a. Up{t,r) = iTr„{^; r^A^''} 

(5.93c) 

Jl = iTrr^vls Sf9. J- = iTr„{s ; CpJ^p], 

(5.93d) 

The velocity operator s is defined in Eq. (5.65) and the 
notation (5.60) is used. 

We now add Eq. (5.93a) to Eq. (5.92) and subtract 
Eq. (5.93b). According to Eq. (5.87) all the collision 
terms as well as the commutators cancel and we obtain 
the energy balance equation [compare with Eq. (2.8)]: 

dtutot + V • il^^ =jE- (5.94a) 

utot{t, r) = u,it, r) + Up{t, r) - Ug{t, r); (5.94b) 
jltM=jl{t,r) +f^(t,r)-jlit,r). (5.94c) 



Ue{t,r) 

we find 



lip Adti - dt^){g{ti,t2,n,r))r 

4 tl — ^t2 — ^t 



ITlVVp 



4 tl^t2^t 



lim {dti - dt2){ng{ti,t2,n,r))n, 



(5.91) 



dtUe + V-jl=j-E 



(5.92) 



+ iU7T lim {dt,-dt,){St{g,J^P,J^^})r 

1 1 — ^t2 — *t 



The expression in the left hand side of Eq. (5.92) has the 
form of the continuity equation for the energy current of 
the electrons and the first term in the right hand side 
is the Joule heat acting as an energy source. The last 
term in the right-hand side indicates that the electron 
system by itself is open, due to the energy exchange with 
the electron-hole pairs. As we discussed in Sec. II A this 
means that in the definition of the conserved energy and 
energy current densities the contribution of these degrees 
of freedom must be taken into account. To accomplish 
this task, we premultiply Eqs. (5.59a) and (5.59b) by 
and £.P respectively. Using Eqs. (5.25), (5.54) and (5.59c) 
and taking the trace Tr„ [see Eq. (5.85)] of both sides, 
we obtain: 



=-Tr„/:9 St" , 



(5.93a) 
(5.93b) 



Equations (5.94), (5.91), (5.93c) and (5.93d) constitute 
the main result of this subsection. They define the con- 
served currents in terms of quantities to be fomid from 
the kinetic equations. It is important to emphasize that 
the conservation laws thus found are exact (at one loop) 
in the sense that no approximation has been made be- 
yond the usual Fermi-liquid theory: in detail, no gradi- 
ent or harmonic expansion has been made and no time 
or space (except those suppressed by a factor of q/pr) 
Poisson brackets have been neglected yet. 

Within the same accuracy with which the kinetic equa- 
tions (5.64) were derived, performing the Wigner trans- 
forms (5.61)-(5.63) of Eq. (5.90), (5.91) and (5.93d) we 
find 

j = evpy j de{nf{e,n;t,r))n (5.95) 
for the electric current density and 
jl=VFi' j dee{nf{€,n;t,r))n (5.96a) 

(5.96b) 

S'"/ -0^{-"^'i^^l)N'{oJ,q;t,r)U (5.96c) 

for the energy current densities, in agreement with 
Eqs. (3.14)-(3.15). 

What remains to be done is to reduce the found ex- 
pressions to the usual form of the quantum Boltzmann 
equation. This is the subject of the next subsection. 
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G. The quantum kinetic equation 



Here and below the short hand notation 



After the Wigner transforms (5.61)-(5.63b), Eq. (5.71) 
becomes 



dt + v-V + ev-E 



f{e,n;t,r) = St{f,NP,N3}. 

(5.97) 

The coUision integral is the sum of the inelastic and elas- 
tic parts 

St{/, NP, NS} = Sti„{/, NP, N^} + Stei{/}. (5.98) 

The elastic part is obtained by adding the transform 
of the "bare" collision integral [the first term on the 

^el 

right hand side of Eq. (5.72)] to the transforms of Sti , 
Eq. (5.75c), and St_, Eq. (5.76b). The inelastic part 
is given by the transform of St^ , Eq. (5.75b), plus the 
transform of St+, Eq. (5.76c). However for the sake of 
compactness we will not make such a distinction between 
elastic and inelastic contributions and, using a notation 
resembling that of section VE, we write the collision in- 
tegral in the following form: 



St{/, NP, N^} = Strf + Sti{/, iV", NS} 

+ St_{/} + St+{/,7V'',iVn, 



(5.99) 



where the first term on the right hand side is the trans- 
form of the "bare" collision integral and the other terms 
are given below. 

Keeping the elastic and inelastic parts of the colli- 
sion integral (5.73) in a single formula, the corresponding 
contribution is obtained by first substituting Eqs. (5.58) 
and then performing the Wigner transforms (5.61)-(5.63) 
[and using their property (5.83)]. We decompose the re- 
sult into distinct contributions due to the two bosonic 
degrees of freedom: 

Sti(e, n; r) = Stj ''(e, n; r) — St^ ^{e,n;t,r). 

(5.100) 

As usual the collisions with the "ghost" particles enter 
with the opposite sign. In terms of the combination 
(5.66) of distribution functions that we denoted by T 
and the vertex (5.28), these contributions are (we sup- 
press the spectator arguments t and r in both sides of 
the equations): 

^e-p 1 f duj 1 f d'^q f dn2dn^dn4, 

(^'"^) = W2^^j/(2^j/ 

X {li2 ['Cg4T4i;2i(e. w> g) + Tg4.2i(e, w, q)Ll^] 

+72l[^34'^42;2l(e>'^.9) + '^34;2l(e>'^.gK42]} 

(5.101a) 

1 fdujl f d'^'q f dn2dn3dni 



Sti (e,ni 



ly J 2ttuj J {2Try 

X {7?2 ['^34T^41;2l(e, t^, Q) + ^§4.21 (f, UJ, q)£-l 
+72l['^^34TC"42;2l(e>'^,9) + '^34;21 (e> 9^42] } 



-41] 



= £"(a;, q; ni,nj), a = g,p 

is used. It is readily seen that this contributions coin- 
cide with the local electron-boson collision integral of 
Eq. (3.19b). 

Proceeding as above, the transform of Eq. (5.76b) is 



St_(e,ni) = Stei(e,ni) -|- St_,;(e,ni) 

d'^q f dn2 



(5.102a) 



- 2„ fdcj 1 
bte; = -Re / 



(27r)° 



d.n(i 2 5 

713746 



X [CP - [/(e - «6) - fie - to, Ui)] 



u J 2n cu^ J i2nr J 



(5.102b) 



drio . . . dn, 



1^5 



-7?37l6 



2/(e, n4) - /(e - uj, 714) - f{e + uj, 714) - (".4 ^ n-e) 



/ 



X / dei 



/(ei,ni)[l - /(ei - w,n3)] + (ni n^) 



X [C - ^-^W- [^'14 - ^'m] (5.102c) 

Equation (5.102b) is the (singlet part of the) elastic 
electron-electron collision integral, Eq. (3.20b). To ob- 
tain Eq. (5.102c) in the given form, we used the analytic 
properties of the propagators and performed the change 
of variable ei ^ ei + lu in some of the terms. 

Finally we perform the transform of Eq. (5.76c) and 
obtain: 

^ ^ e-e ^ ^ 

St+(e, ni) = St„ (e, ni) + St+.„(e, ni) + St+,i(e, ni) 

(5.103a) 

The first term is given by (the singlet part of) Eq. (3.20d). 
The second term is 



^ e-p ^ e-g 

St+^„ = St_|_ „ — St^ „ 



(5.103b) 



with St_|_„ given by Eq. (3.19d) excluding the last line. 
The third term is: 



d'^q 



dn2 ■ ■ ■ driQ , - 

713746 



(5.101b) 



2uJ 2nujJ {2-kY J 

2/(e, m) - f{e - ui, n4) - /(e + w, m) - {m ne) 

X [C {2NP + 1) - (27V9 + 1)]52 [£§4 - ^14] (5.103c) 

Adding Eqs. (5.103c) and (5.102c) we recover the last 
line of Eq. (3.19d) and the (the singlet part of the) local 
electron-electron collision integral, Eq. (3.20f). This con- 
cludes the derivation of the quantum kinetic equation for 
the singlet channel. 



H. The triplet channel 

The inclusion of the interaction in the triplet channel 
is straightforward; in the Eilenberger equation (5.7) we 
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add to the left hand side the term: 



(5.104) 



where are the Pauli matrices [i = x, y, z] and the fluc- 
tuating field is a 3 component vector in the L = 1 an- 
gular momentum space. Therefore all the triplet channel 
propagators, the polarization operators and the density 
matrices should be treated as 3 x 3 matrices; for example 
we have 

[DoL (1,2) = d{tu)Oij 

■' V 

[cf. Eq. (5.15)] and the retarded polarization operator is 
given by 



n^(l,2) = i.|^i2 + Jtr 



<5ff (ii,ii,ni,ri) 



(5(?i+j(t2,r2,n2) 



where the trace is over spin indices. 

The trace of triplet channel operators includes the sum 
over the indices In the absence of the magnetic field, 
all the operators are diagonal, e.g.: 



and the trace results in extra factors of 3 in compari- 
son to the singlet channel. The derivation can therefore 
be repeated with simple modifications and it gives the 
quantum kinetic equation presented in Sec. III. We note 
only one main difference in the derivation for the triplet 
channel: the gauge transformation, which has the form: 

does not commute with the interaction term (5.104). Ad- 
ditional second order terms arise due to commutators of 
the Pauli matrices; however these terms vanish in the one 
loop approximation and we can neglect them. 

In the next section we will use the quantum kinetic 
equation to calculate the interaction corrections to trans- 
port coefficients and specific heat. 



VI. DERIVATION OF TRANSPORT 
COEFFICIENTS AND SPECIFIC HEAT 

In this section the calculation of the transport coef- 
ficients for quasi one-dimensional, two-dimensional and 
three-dimensional systems is presented; the evaluation of 
the interaction correction to the specific heat is in the 
final subsection. 

In order to calculate the currents in the presence of an 
external field (electric field E or temperature gradient 
VT), we need to solve the kinetic equations. We assume 
that the external fields are weak, i.e. 

eELT < T, VTLt < T 



and the temperature length was defined in Eq. (5.53). 
These conditions ensure that the deviations from the 
equilibrium distribution functions are small and we can 
solve the equations by iteration. 

In the lowest order, the distribution functions should 
null the collision integrals in Eqs. (3.4) and (3.6), and the 
sought corrections (5/, 6N'^ are linear in the electric field 
or in the gradients of the distribution functions. In other 
words, we look for a solution to those equations of the 
form: 

/(e,n;r) = /;^(e;r) + 5/(e,n;r); (6.1) 
7V"(u;, g; r) = Np{uj; r)i + SN^^lu, q; r) ; 

where the Fermi and Plank distribution functions (2.2) 
depend on the spatial coordinate only through the tem- 
perature T(r). For compactness, we will consider explic- 
itly the singlet channel only and we will indicate how to 
include the triplet channel. 

Let lis start from the electronic part of the kinetic equa- 
tion. The bare impurity collision part St,- is larger than 
the other terms, thus it suffices to calculate the latter in 
the first order of perturbation theory. Considering short 
range impurities, so that r is independent of the scatter- 
ing angle, we find: 



5f = 5fo + 6h; 
Sfo = TV- \eE- 



(6.2) 



eVT 



dfF{e) 



T I \ de 

6fi = rSSt {fp + <5/o, Np + SNP, Np + SN^} 

where SSt is the linearized collision integral. Note that 
according to the discussion of the conservation laws in 
Sec. Ill, we need to consider only the local electron- 
boson contribution, Eq. (3.19b), and the elastic electron- 
electron one, Eq. (3.20b). 

Expression (6.2) is to be substituted in Eqs. (3.14) and 
(3.15b) to find the electric current and the electron com- 
ponent of the energy current. Integration of the Sfo term 
is straightforward. Due to the structure of the collision 
integrals (3.19)-(3.20), the integration over e can be per- 
formed before the lu and q integrations in the Sfi term. 
For the combination of distribution functions entering 
into Eq. (3.16a), we find 

(5Ty,fei(e,a;) = 6N{u!, q;ni,nj) {/^(e) - /^(e- w)} 
+ fld5{n^j)Np{u>) {Sf{€, Uk) - 5f{e - w, nj,)} 

+ Q.d5{ninj) (5/(e,ni)(l - /^(e - w)) 



which with the help of the identities 

.dfpie-oj) d 



defpie)- 



de 



dui 



ujNp{uj) 



de 



(6.3) 



(6.4a) 
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I 



3 5w 



dee(e-a;)/F(e) 



r£~A 

3 dtjj 



LuNp{uj) 

de 

(jjNp{ijj) 



1 ^ dNp{(jj) _ 
3 doj ' 



(6.4b) 



1 ^dNpiu))^ 
6 9a; ' 



yields 

^ y de[(5Tij,fc;(e,a;, q) - ^Tij,fe;(e, -w, -q)] 

^ y rfe e [(5Tjj,fe; (e, w, q) - STij^kl (e, -t^, 
= — ^5N{lo, q; Ui, Uj) — (n, + nj.) w 



(6.5a) 



(6.5b) 



vftVT 
+ — Tf; {ni-nk) 



du) 



12 Olo 



where 



Nf 



{rii, rij) = Np{u))Cld6{ninj) . (6.6) 



Here we retained only the odd in uj contribution because 
the even part vanishes after the u) integration in the rel- 
evant collision integral, see Eq. (3.19b). 

The combination of the distribution functions entering 
the elastic collision part, Eq. (3.20b), gives 

i J defpie) [6f{e - n) - df{e + to, n)] (6.7a) 

F) r 1 

LoNp{uj) . 



-evprE ■ n 



duj 



e/F(e) [Sf{€ -uj,n)- Sf{e + uj, n)] 



(6.7b) 



vpTn ■ VT 



du! 



6 du! 



where once again we retained only the odd in u part, 
non- vanishing after the oj integration in Eq. (3.20b). 

Using Eq. (3.14), we find the electric current j = aE 
with the conductivity tensor & given by: 



diij 



OU) 



(6.8) 

where Oo = TVpe'^v/d is the Drudc conductivity. In- 
dicating spatial indices with jj,!/ = 1, . . . ,d, the elastic 



kernels 5*;'^ and 5^^, which originate from Eq. (3.19b) 
and Eq. (3.20b) respectively, are: 

11 d r d'^q f dnidn2 

'^-^") = ;^ 7(2^ 710^"^""^^ 

Re - £fi] + Re - C'^^] - 2Rc [Cl^ - Cf,] 
ci2/ N 2rf f d'^q fdnidn2 . p g, 

dr f d'^q fdni...dn6 Q 5 

= y (2^ y ^f]^^-^- 

{nl^,n3^ - nz^ng^) Re| [C^q - rfg] [£^3 + Cf^] | 

(6.9) 

where for compactness we kept the singlet channel cor- 
rection only; inclusion of the triplet channel contribution 
is straightforward'*'^. We show in appendix F that our 
expression for the conductivity coincides with the one of 
Ref. 17. It is natural that the conductivity does not in- 
volve any bosonic distribution function [cf. Eq. (6.5a)], 
because the inelastic electron collision with such bosons 
changes the energy of the electron but not the direction 
of its motion. 

In contrast, even the electron contribution to the ther- 
mal conductivity tensor k, such that jj^j = —k^T, is 
sensitive to the bosonic distribution functions. We rep- 
resent the total thermal conductivity as 



+ 6k + k'' 



(6.10) 



The first term in this expression obeys Wiedemann- Franz 
law even with the interaction to the conductivity in- 
cluded, i.e. kyyp = haT with the Lorentz number given 
by Eq. (1.1). The second term represent the (electronic) 
correction to the Wiedemann- Franz due to the energy de- 
pendence of the clastic scattering and due to the inelastic 
electron scattering on the bosons. Finally, the third and 
the fourth term represent the contribution of the p and g 
bosons to the thermal transport. These additional con- 
tributions are given by: 



6k = Sk, 



el 



Sk,, 



(6.11a) 



12 dw 



(6.11b) 



+ Vf 



d'^q 



^3 dNp 



fdu) f 



4 du 
dnidn2dn3 



{2nY J {^df 
{ Re [C{^5,NP^ + CP^^SuNP, - C{^5,N^^ 



-Re 



r9 K pj9 _i_ r9 X pj9 _ r9 X m9]\. 

(6.11c) 
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dui f d'^q fdnidn2 r.^, x ^^a 1 



2Tr J {2TTy J (Qd) 



with 



(6.12) 



S^N^j = [dN"{uj, q; m, n,-)] • (6.13) 



Equations (6.8)-(6.12) arc the complete expressions for 
the electric and thermal transport coefficients. To obtain 
the explicit result one has to solve Eqs. (3.6) to find the 
distribution functions dN". We will do so by restricting 
ourselves to the diffusive Tr <C 1 regime, except for two- 
dimensional systems for which we will consider the ar- 
bitrary temperature range'*^. Moreover for the Coulomb 
interaction we will consider the unitary limit (for infrared 
finite momentum integrals), which will enable us to drop 
all the terms that depend ondqF. 

Prom now on we will retain only the zeroth harmonic of 
the Fermi liquid constants, which wc denote symbolically 
with Fq . For the singlet channel this means 



and for the triplet channel 



(6.14a) 



We remind that — and that for the Coulomb inter- 
action the potential is: 



V{q) 



2_ 



eMn 



{qaf 



d = 3 

d^2 
,d=l 



(6.14c) 



where a is a length of the order of the quasi one- 
dimensional wire width. 



For the distribution function TV, the above is only the 
starting point for the iterative solution: 



6NP =SN° + 6N^. 
The equation for 6N^ is: 



(6.17) 



+ i 



(jO dNp{uj) 
T duj 
pp 

'l + pp' 



d FP 
dql + FP 



2{st^;57V^} . 

(6.18) 



In the diffusive limit Tt <C 1 the (first iteration) solution 
would be of the form 



du) 



for a vector V with magnitude of order one. However 

contributions from frequency lo larger than temperature 
T are exponentially suppressed, i.e. ujt ^ Tr <C 1; there- 
fore 6N^ can be neglected in comparison to SN^ . Thus 
in the diffusive limit 



T du 



(6.19) 



For the propagators £" the diffusive approximation 
amounts to the substitution 

which leads to the following expression: 



C"{ni,n2) =T{6{rirfi2) - l) + 

-I- (m + n2)u,L^f, + ni^n2i.L^f,^, 



(6.20) 



where the functions Lf depend on w, q only and are ex- 
plicitly: 



A. Diffusive regime 

We consider first the distribution function N^; substi- 
tuting expression (6.1) into Eq. (3.6) we obtain at linear 
order in VT 



-I- Vp 



LO dNp{u)) 
T doj 



+ i 



vq;5NS 



:{st^;5Afs| 



VT • (0^(5 (nin2) ni - ni - ^2) 



(6.15) 



with Np defined in Eq. (6.6). The (exact at this order) 
solution for 5N^{w,ni,n2) is: 



5NS = SN'' = vpTUi ■ VT 



T du) 



(6.16) 



L^^,=-dTDq^q,L^. 



(6.21a) 

(6.21b) 
(6.21c) 



Here D = rvp /d is the diffusion constant. These formu- 
las are valid for u), Dq^ <^ 1/r. 

Within this approximation, we have 

This means that in both a and 5k sec Eqs. (6.8) 
and (6.11) - we can neglect the contributions of the S 
kernels [note that by inserting the solution (6.19) into 
Eq. (6.11c), Skin is given by twice the first line of that 
equation] : indeed the leading contribution is given by the 
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kernel £ in Eq. (6.11b). This kernel has the approximate 
form 



i 1 



nvT duj J {2ttY 



Re 



Lm,^-L'2,.) ■ (6.22) 



Finally the bosonic contributions [cf. Eq. (6.12)] to the 
thermal conductivity (6.10) can be written as^^ 



dNp 



4 Olj 



with 



(27r 



(6.23) 



(6.24) 



The next step is to evaluate the momentum integrals; 
first we give the results for the short range interaction de- 
scribed by the (zcroth harmonic of the) Fcrmi-liquid con- 
stant Fo and then we indicate the modifications needed 
to account for the long range part of the Coulomb inter- 
action in the singlet channel. The triplet channel contri- 
butions are obtained by multiplying the found results by 
three and identifying Fq with Fq . 

For the elastic kernel we find: 



f/i!/(w) 



d CFo {2ttY 

1 



1 



cos 2^ 



|-±(l + F„-(I + F„,M) 

(6.25) 



The expression for the Coulomb interaction is obtained 
by taking the unitary limit ^ +oo. 
The result for the bosonic kernel is: 



2 an {2n)'i 
1 



1 

u 



cos 



l-(l + Fo)'-3 



(6.26) 



In d = 3 the limit Fq +oo gives the correct formula 
for the long range contribution, but for d = 1 the limit 

diverges. By retaining the full form of the interaction 
potential however such infrared divergence is cut off at 
the inverse of the screening radius. With logarithmic 
accuracy the result for the Coulomb interaction is found 
by substituting [. . .] —ak\n^^^{Dk^/\u}\), where a is a 
length of the order of the wire width and fc^ = Ane^v is 
the square of the inverse screening radius [in the bulk]. 

We can now proceed with the calculation of the inte- 
grals over oj in Eqs. (6.8), (6.11) and (6.23) using the 
identity: 

j dij u'^d^Npiuj) = -2r'"C(m)r(m + 1) (6.27) 



Here (^{x) is the Ricmann zcta function, whose values at 
the points relevant for our discussion are: 



1 

' 2 ' 



C(-i)« -0.208 , C(0) 
C(i)« -1.460 , C(i)« 2.612, 

C(2) = ^« 1.645, C (I) « 1-341, 
and T{x) is the Euler gamma function, with values 

r(i) = 0F , r(i) = i, 
r(|) = i,/^, r(|) = f^. 



r(3) = 2 
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Performing the final w integrations we obtain: 



cr =(Td + e 



{2ttY d 
d-A 1 
2 - cos ^ 



(6.28) 



5k = 



cos 



Tzd 



d 1 / 1 ti 

---(l+Fo-{l+Fo)'-^ 



(6.29) 



1 



DC 



4:{2nY \ D J ' V2 
1 



1 ir( - + 2 



cosf 



l-{l+Fof-i 



(6.30) 



where in the absence of the magnetic field tr^jy = aSfii, 
and the similar relation holds for the thermal conduc- 
tivity. According to the previous discussion, for the 
Coulomb interaction the correct expressions are given 
by the limit Fq — > -|-oo, with the exception of the 
term k'' — in d = 1 for which, with logarithmic ac- 
curacy, the result is obtained by substituting [. . .] — > 
— afcln^^^(D/c^/T). The final answers for the corrections 
to the thermal conductivity are given in Eqs. (4.2)-(4.3). 

We note that for d = 3 the oj integration in a and 
KwF is actually ultraviolet divergent. This divergence 
can be incorporated as a renormalization into the Drude 
results; this renormalization however does not invalidate 
the Wiedemann-Franz law. 



B. Two-dimensional system 

In order to evaluate the interaction corrections for the 
whole temperature range, we need the exact form of the 
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propagators. In two dimensions they are: 

q; ni, 712) = {n^ri2) £o(w, q; rii) 



(- 



(6.31) 



where 



jCo{w,q;n) = — — — — , 

—iLO + IV ■ q + l/r 

C=^{-iLO+l/Tf + {vFqf. 



(6.32) 



Note that now the variables u) and vpq are hmited only 
by the Fermi energy Ep- 

As before, we need to find the non-equilibrium correc- 
tions SNP'^ to the bosonic distribution fimctions. These 
are given again by Eqs. (6.16)-(6.17), but to obtain the 
thermal conductivity on the arbitrary temperature range, 
we calculate exactly (at linear order in VT) the solution 
to Eq. (6.18): 



6N^ = N 



qVT 



vpq" 



ivpq ■ (rii - 722) + 2/t 



X{ni)ni — X{n2)n2 



(6.33) 



where the bar indicates complex conjugation and we in- 
troduced the following quantities: 

- UJ^dNp Fq ^ {nq)q 

N = vpT—— — z ^ , n = n ■ 



X{n) 



T du) 1 + Fo 



«('^) = \ {'"FqY + ( -iq v 



We will not calculate the corrections to the electri- 
cal conductivity, which would reproduce the results of 
Ref. 17, as shown in Appendix F [see also Ref. 18 for 
the generalization to the arbitrary disorder]. For conve- 
nience in the calculations, we separate in 5 kin and k'' the 
contributions due to ^A''^ and 5N^ [cf. Eqs. (6.11c) and 
(6.23)]: 



5kin =6k°„ + 6kl„ 



Lj^ dNp 



4 duj 



VF 



cko 
2^ 



cPq f dnidn2dnz 



(27r) 



{^2 



(6.35a) 

(6.35b) 

(6.35c) 



dNf 



4 duj 
dnidn2 



(6.36a) 
(6.36b) 

(6.36c) 



where, as in Eq. (6.13), we indicate with d^N^ the varia- 
tional derivative with respect to the temperature gradient 
and 



d^q f dni 



-ni^ni^Re [C'{-y - £f J . 



nujv J (27r)2 J ^2 

(6.36d) 

Expressions (6.8) for the electrical conductivity and 
(6.11b) for the elastic correction to the thermal conduc- 
tivity and the definition (6.9) for the kernels remain un- 
changed. The momentum and angular integrals in the 
latter can be exactly calculated and for the singlet chan- 
nel in the unitary limit we find: 



e^V 1 



2 — 2u}tH{u)t) arctan 



1 



2 ^('^-)jin^+(^ 



(6.37a) 



-{ujTfH{ujT)ln2 



-^r-sgnw 



o-D27r2 



-|--sgna; -|- u)tH{u)t) 



(6.37b) 

I (2H{ujt) - Ij arctan J- (6.37c) 

-hln2 



1 / 1 



with the function H defined as: 



H{x) = 



(6.37d) 



To perform the momentum integral in Eq. (6.36d) we 
must keep the full form of the propagator in order to 
avoid the infrared divergence that we would obtain in 
the unitary limit: 



BU^) = - 



(Td 277^ (jO 



-< In 



-H TT W T 



- In [1 + (u^'t)^] — WT arctan | 



(6.37e) 



where k = is the inverse screening radius. Next, we 

calculate the angular integrals in Eqs. (6.35c) and (6.36c) 
as well as the angular part of the momentum integrals 
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and obtain: 



d9q f dnidn2dn^ 
27 



3 

1 1 



~ 4.^^'''''\ {vqY CC 



Cb' Cb' 
+ 



+ 



{vqy 



C-b' 



C-b 
C-b'' 



C-b 
1 

CC 



[C-b'-C + b' 
C C 



C-b 



C-b\ 
(6.38a) 



The final integration can now be performed within the 
logarithmic accuracy; we find [cf. Eq. (4.4b)]: 

A.. = -^r(Tr)^ln(^ 



+ 



where 



(6.41) 



c 



C-b 



1 



Cb' Cb' 

+ 



C-b^ 
C-b' 



C-b C-b 



1 1 

c + c 
c- 



(6.38b) 



b' 



C 



The function C was defined in Eq. (6.32) and b' is given 
by the limit Fq — > +oo of b in Eq. (6.34), where N is 
also defined. The remaining integrals over the magnitude 
of the momentum can be evaluated approximately; the 
result can be written as: 



dNp 



4 duj 



(6.39a) 



with 



Ef 



■ arctan wr — ujt In 



(6.39b) 



In the above kernel the first term in curly bracket orig- 
inates from only - as discussed in Sec. VB no long 
range terms can be present in the electronic contribution 
to the thermal conductivity. Note that the second term in 
the above expression is beyond the logarithmic accuracy 
of our approximate calculation and must be dropped. 
Similarly, most of the terms in the other kernels can be 
neglected and collecting the logarithmic contributions: 



1 



27r2T 

2 



dio 



dNp 



(a;r)^ln(:^) + 



7 



1 + {u^Tf 

16 1 



In 



^(c.r)2-g-ry^^^ 



ln( H 



1 



(6.40) 



gi{x)=— \ - 



2^'[- 

X 



3i{x)=-^x 



3^1(0.) -3, 



(6.42a) 
(6.42b) 



and tp' is the derivative of the digamma function. Since 
the asymptotic behaviour of gi (x) is 



.x < 1 

x > 1 



(6.42c) 



both these functions tend to 1 as Tt 0; therefore in the 
diffusive limit the main contribution is T\si(Dk^ /T)/V2. 
On the other hand, for Tt 1 the first term in Eq. (6.41) 
is the dominant one. 

Turning to the triplet channel, for simplicity we will 
restrict ourselves to the limiting diffusive and quasibal- 
listic cases, although one can extend the calculation to 
the whole temperature range, as done in Ref. 17 for the 
electrical conductivity. 

In the diffusive limit Tt <C 1 we know from our previ- 
ous analysis that we can discard the B"^ term as well as 
the S terms. The relevant kernels are then'''^: 



^/ni.(w) = - 



CFd 



1 



1 



ln(l + Fo") 



ln(l + Fo-) 



(6.43) 



(6.44) 



which substituted into Eqs. (6.11)-(6.36b) give 



T 
18 



1 



1 



ln(l + Fo- 



T 
12 



(6.45) 

In the opposite limit Tt ^ 1, the main contribution 
comes, as for the singlet channel, from the logarithmic 
divergence at large momentum in the kernel B^: 



"(7D27r2 



-In 



Ep 



(6.46) 

Then, the correction to the thermal conductivity is: 



where we defined the singlet and triplet channel correc- 
tions as: 



A.. = -^r(Tr)2ln(f: 



per 



l + i^r? 



(6.47) 



[k - K„p]^^ = {Aks + 3AKt) 6^,^ 



which concludes the derivation of Eq. (4.4c). 
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This correction to the thermal conductivity (and the 
corresponding one in the singlet channel) is the inelastic 
processes contribution to the energy relaxation rate^^. 
In the clean system, such inelastic processes cannot relax 
momentum (because of Galilean invariance) and hence 
they do not affect the electrical conductivity, but they 
can contribute to the energy relaxation rate T^. In the 
kinetic theory, the thermal conductivity can be written, 
up to a numerical coefficient, as: 

K oc TEf/T^ 

and the rate is given by the sum of the rates for the 
relevant processes, namely the electron-impurity and 
electron-electron scattering rates: 



with Timp = 1/r and 



T 



Fe-e = a— In I — 



Ep 



wr 



ith^e 



)L (6.50) 



As before;, we will consider explicitly the singlet channel, 
short range interaction in the zeroth harmonic approx- 
imation for the Fermi liquid constant, which we denote 
with Ff). The results for the long range interaction in 
the unitary limit are obtained by letting Fq +co. For 
the triplet channel one must substitute Fq with Fq and 
multiply by an overall factor of 3. The final answer with 
the correct coefficients is given in Sec. IV B. 

In the diffusive limit, to which we restrict our attention 
for d= 1,3, we have: 



1 



(6.51) 

and the functions Lq were defined in Eqs. (6.20) and 
(6.21). After the integration over momentum we find: 



Here a is a constant whose exact value is irrelevant for 
our argument. In the limit {T'^/Ep) \n{EF/T) < l/r we 
can expand the expression for the total rate, substitute 
the result into the above formula for k and obtain: 



Ep 



K a TtEf - aT{TTY In ( ^ 



The first term on the right hand side is the usual Drude 
result for the thermal conductivity and the second term 
has the form of the correction (6.47). Note that in the 
opposite limit (clean system) one would find: 



in agreement with the result of Ref. 29. 

C. Specific heat 

Recalling our discussion on the structure of the kinetic 
equation in Section II, we write the total specific heat as 
the sum of the electronic and bosonic contributions: 



Cy — Cy -\- SCy^ 

Q dUe TT^ 
^Cy = -g^(up-Ug^, 



(6.48a) 

(6.48b) 

(6.48c) 



where, according to Eq. (5.93c), the bosonic energy den- 
sities [in the equilibrium (3.1)] are 



j dujujb"{ij)Np{w) (6.49) 



b''{uj)-b%uj) = 



cos f (d - 
sinf (d- 2) 



1 

4w 



2) 



(6.52) 



which inserted into Eq. (6.48c) gives: 



SCy 4 



2 (27r) 

cosf (d- 2) 
sin I (d- 2) 



1- 



(6.53) 



The relevant numerical values for the zeta and gamma 
fmictions are given after Eq. (6.27), which has been used 
to perform the oj integral. 

For = 2 we can keep the full form of the propagators 
to find the singular contribution to the specific heat at 
arbitrary value of Tr: 



Re 



27r J (27r)2 
i-iLU + l/r) / 1 



^0 



1 



F^C 



C-b C-l/i 



(6.54) 



with C and b defined respectively in Eqs. (6.32) and 
(6.34). The first term in the integral is formally diver- 
gent at |g| DO and this divergence gives a linear in T 
contribution to the specific heat which does not depend 
on disorder. This term must be disregarded as all the 
linear terms are included in the definition of the effective 
electron mass - taking it into account here would be a 
double counting. To regularize the integral, we replace 
— ^ in the first line. 
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Performing the momentum integral we obtain: 



^0 



-t\^\ 



(6.55) 



The final answer for the correction to the specific heat is 
then: 



5Cv = - 



1 T 1 
T2Z)TTJ^ 
1 



ln(l + PP) 



47r2 



(27-3)C(2)-2C'(2) 



^0 



(6.56) 



where C(2) « 1.645, ('{2) « 
In the quasi ballistic limit r 
relevant: 



Sc. = --C(3) 



-0.938 and C(3) w 1.202. 
+00 only the last line is 



pP 
^0 



This correction to the specific heat has the same form 
found for two-dimensional Fermi liquids^*^. 

As discussed before, the long-range interaction can be 
accounted for by taking the limit Fq — > +00, while the 
triplet channel contribution is three times larger - see 
also Sec. IV B. 



Rcf. 11 for the diffusive regime and in Ref. 17 in the 
ballistic and crossover regime. 

The neutral bosons, however, are crucial for the ther- 
mal properties of the system. Namely, their contribu- 
tions to the energy density are responsible for the non- 
analytic corrections to the specific heat, see Eqs. (4.6) 
and (4.8). Our kinetic equation approach reproduces 
the results for the interaction corrections to the spe- 
cific heat previously calculated within the equilibrium 
diagrammatic technique^^ . Moreover the neutral bosons 
contributions to the energy current cause the violation 
of the Wicdcmann-Franz law, sec Eq. (4.1) and the dis- 
cussion that follows it. The violation is stronger for low 
dimensionality systems {d = 1,2) in the diffusive regime, 
sec Eqs. (4.2) and (4.4b). Other effects contributing to 
the violation of the Wicdcmann-Franz law are the energy 
dependence of the electron elastic scattering and the in- 
elastic scattering of the electrons on the neutral bosons. 
The latter effect was found to be relevant in the quasi- 
ballistic regime Tt h for two-dimensional systems, see 
Eqs. (4.4). 

The violation of the Wiedemann-Franz law was inves- 
tigated before in the diffusive regime in Rcfs. 20 and 23 
within the "quantum kinetic equation" approach and by 
Kubo formula in Ref. 22. Ironically, even though the 
forms of the energy current operator used in those refer- 
ences are wrong, the final results for the thermal conduc- 
tivity are consistent with our Eqs. (4.2)-(4.4). We think 
that this agreement is accidental. 
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VII. CONCLUSIONS 

Locality on the scale determined by the temperature 
and the validity of the conservation laws are two main 
requirements for a proper kinetic description of any sys- 
tem. In the present paper we derived such a description 
for the interaction effects in disordered metals [assuming 
that the clean counterpart of the system is a stable Fermi 
liquid] . 

We showed that this description requires the introduc- 
tion, along with the usual fcrmionic quasiparticlc distri- 
bution function, of additional bosonic distribution func- 
tions. These neutral bosons are of two types: (i) the 
ones describing the oscillation in charge density (singlet) 
or spin density (triplet) and (ii) fictitious (ghost) bosons 
which prevent over-counting the degrees of freedoms 
(electron-hole pairs) already included in the fermionic 
part. The obtained conservation laws together with 
gauge invariancc allowed for the unambiguous definition 
of the corresponding electric and energy currents. 

For the electric transport the neutral bosons are not 
important and our description reproduces the known re- 
sults for the correction to the conductivity obtained in 
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APPENDIX A: ON THE CORRECTION TO THE 
THERMODYNAMIC POTENTIAL 

Standard analytic continuation of Eq. (2.17) gives: 

=/sr°">(^)-'^>»('+7ft") (^^) 

In the second line we used the operator notation see 
Eqs. (3.3) and (3.13) - which gives the correct generaliza- 
tion for the momentum dependent Fermi liquid parame- 
ter. Substituting the transform of the explicit expression 
(5.43a) for the polarization operator, the argument of the 
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logarithm can be rewritten as: 



(i+f) {i^y +- 



+ F 



According to the definition (3.11), the term in square 



brackets is y^'^j ■ Using the property 

Tr In (^AB^ = Tr Ini + Tr InS 

and since In ^1 + does not contribute to the imagi- 
nary part, we conclude that 



ImTr In 1 1 + j = -ImTr 



In/;" -In £9 



Substituting this identity into Eq. (Al) we finally obtain 
Eq. (2.18). 



APPENDIX B: ON THE MICROSCOPIC FORM 
OF THE ENERGY CURRENT OPERATOR 

The action entering into the partition function that 
describes the electron gas in the presence of an external 
electric field is 



dtd'^r 



(Bl) 



with the condition: V x A = 0, which ensures the absence 
of the magnetic field. The variables {t, r) on which all the 
fields depend have been suppressed. The gauge invariant 
part of the Hamiltonian for the non-interacting system 
has the usual form: 

V'^4iV'=^(^V + A)^t(_iv + A)^ (B2) 

where Vimp is the impurity potential and the potentials 
(fi and A describe the external electric field: 

eEext = -Vif + dtA 

As usual, the charge conservation law 

dtp + V-j = (B3) 

follows from the requirement of gauge invariance, with 
the charge and current densities given by: 



3 = 7^ [V"^ (-«^ + A)ilj + (iV + A) VV] 

^TTt 



(B4) 



The invariance of the action upon the replacement 

-ip^t, r) ip{t + a{t, r),r) 



[and the similar replacement for t/j^] is the basis for the 
derivation of the energy conservation law. A straightfor- 
ward calculation gives 
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= = d^{i>^Hgi^)+^dt{^^^) 



2m 



V • [dti^\-iV + A)V - (iV + A)V^9t^] 



(B5) 



where the prime means that the derivative acts on tp, ip'^ 
only. By adding and subtracting terms proportional to 
in the last bracket and to dtA in the first term, we find 
the energy conservation law: 



dtUo + V-j°=j- Eea:t - [dtp + V • j] 



where 



U0=1p^Hgi1p 

-{iV + A)^\idt 



ip)xpH-i'V + A)'ip 
ip)tp 



(B6) 

(B7a) 

(B7b) 



The first term on the right hand side of Eq. (B6) is the 
usual Joule heat; the last term in that equation is not 
gauge invariant, but it vanishes because it is proportional 
to the continuity equation (B3). 

We now consider the generalization to the interacting 
case. The Hamiltonian now contain an additional term: 

d:^ritp^tp{r)V{r -ri)'p'^'p{ri) 

where V{r) — e^/|r| describes the density-density 
Coulomb interaction, which can be decoupled by the 
Hubbard-Stratonovich transformation. This amount to 
the introduction of the quantum fields <j) and A. in the 
action by adding the term 



.^](j)^ + \E)i, eEfi = -V(j) + dtA 



(B8) 



and redefining the vector potential as the sum of the 
external and fluctuating ones: 



+ A 



(B9) 



The variation of the action with respect to the fluctuating 

potentials results in the first and fourth Maxwell equa- 
tions relating the fiuctuating electric field to the charge 
and current densities: 



V-Efi=p, = j + dtEfi 



(BIO) 



where the electric current is defined in Eq. (B4), but the 
substitution (B9) must be performed. 

To obtain the energy conservation law we must con- 
sider the further transformation 

^{t,r)-^ (P{t + a{t,r),r) 
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and the similar one for A~ Proceeding as before we find 
the conservation law: 



du^ „ .? 
1- V • 7 

dt 



Uo 



where uq, j° are defined in Eq. (B7) [with the substi- 
tution (B9)]. Given the form of Eq. (Bll), one might 
be tempted to call u' and the energy and energy cur- 
rent densities so that the conservation law takes exactly 
the same form as in the non-interacting case. However, 
such a redefinition would result in gauge dependent ex- 
pressions for the densities, since the terms the square 
brackets taken separately are not gauge invariant. Hence 
this naive redefinition of the conserved quantities is un- 
physical, since any physical perturbation can be coupled 
only to gauge invariant quantities. To find the gauge in- 
variant definitions, we rewrite the contribution of those 
non-gauge invariant terms as: 
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(Bll) 
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d_ 
dt 



OA 
'dt 



E 



fi 



<f>3 



dt 



E 



fi 



= ^ct>[dtP + V-j] + ^^[p-V-Efi] 

+ ~ {dtA~ V<l>)-Efi + ^ [V(/. • j + dtAdtEfi] 

Here the first line vanishes because of charge conserva- 
tion, Eq. (B3), and because of the first of the Maxwell 
equations (BIO). As for the second line, we use the sec- 
ond Maxwell equation to eliminate the current; in the re- 
sult we substitute the definition for the fluctuating field 
given in Eq. (B8) and we obtain that the second line of 
Eq. (Bll) equals dtE^^i. This enables us to conclude that 
the correct, gauge invariant expressions for the energy 
and energy current densities are: 



-i,^H,ii,+ \El 



(B12) 

(B13) 



. _ 1 

-{iV + A)^p\idt-ip)^ 
where the potentials are the total ones: 

A = Aext+ A, <^ = <^ext + 

We note that these expressions are gauge invariant with 
respect to gauge transformations of both the external 
and fluctuating potentials. We believe, that only such 

quantities can be coupled to the "gravitational field" in 
the Luttinger scheme for the calculation of the thermal 
conductivity^. 

The same final answer is obtained if the interaction is 
decoupled in the "gauge-fixed" form A = 0. In this case, 




(a) 



(b) 



FIG. 4: (a) the energy current vertex for the non-interacting 
case; (b) the additional vertex induced by the interaction. 
The solid lines with arrows arc electron Green's functions, 
the wavy line is the interaction propagator, the dashed lines 
are the "standard" (non-interacting case) energy and electric 
current operators defined respectively in Eqs. (B7b) and (B4). 



which is the most widely used in the literature, there 
are two contributions to the energy current vertex in the 
diagrammatic approach, see Fig. 4: beside the usual ver- 
tex of the non- interacting case, which arises from the 
terms dr'tp^'Vip, there is an additional vertex from the 
(jitp^^tp terms. These vertices were not taken into ac- 
count in Refs. 20,22,23. However, analogous vertices 
were previuosly considered in the calculations of the ther- 
moelectric coefficient with the inclusion of the electron- 
electron interaction in the particle-hole channel^*^'^^ and 
in the Cooper channel^^ and for the electron-phonon 
interaction^^. 



APPENDIX C: ALTERNATIVE 
PARAMETRIZATION 

The operator Hg.h defined in Eq. (5.54) is clearly not 

a standard Hamiltonian. However, one can introduce a 
different definition for the propagator 



— +^^^e-h(-^Vl)-St, 



(CI) 



such that the (new) He-h operator is indeed a Hamilto- 
nian: 

He.h{q) = Pv-qP, St°"* = PStrP, 

1/2 (C2) 



:F=[i + Fj 



[the action of the operator F was defined in Eq. (5.49)]. 
Proceeding as in Section VD, one obtains the following 
expressions for the /C propagators: 



{dtY 



vt^ = - CP - TLPp {dt) 



K 



(C3a) 



{dt)-^ (C3b) 
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The "kinetic equation" for , Eq. (5.64a), remains un- 
changed, while M'^ now obeys: 



dt + = St" {M<',g''] 



(C4) 



with 



St'' {A/"", 5^} EE 2|st™*; A/""} + 2TMp (C5) 
or, after the Wigner transforms (5.61)-(5.63): 



dtNP + {a, VTV"} + i [iJe.h((z), NP 



-H 



St'^{7V^/} (a;,ni,n2) 

^16734 '^52;43(e' '^)-^45 + ^56734 '^15;34(e, w)^42 



= St^{iV'',/} 
(C6) 

(in3 . . . driQ 

(C7) 



Q4 



with the definitions (5.27) and (5.66) for 7 and T. 
We can then proceed as in Sec. VF and obtain the 

conservation laws (5.89) and (5.94); the only formal dif- 
ference is in the definition of the bosonic energy density, 
which now is: 



Ua{t,r) 



1. 



-Tr, 



In the alternative parametrization the formalism can 
be developed with no more difficulties than in the original 
one. However the evaluation of the thermal conductivity 
become cumbersome. Also in the original parametriza- 
tion it is easier to include (at least perturbatively) the 
effcicts due to higher harmonics of the Fermi liquid pa- 
rameters. 



APPENDIX D: DERIVATION OF THE 
ELECTRONS COLLISION INTEGRAL 

The calculation of the matrix collision integral is sim- 
plified by the introduction of two functions A{t, r, n, n) 
and B{t,r,n,n) such that 

^iK{t,n,r)^-iK{t,n,r) ^ ^ a (Dl) 

We remind that K = K+Ik + ^-O'x 



'K 



f 

1 



The collision integral - right hand side of Eq. (5.22) 
in the matrix notation is then: 



1 

27' 



g{n)o{AiK + Bal)^J{h){AiK + Bal] 



where the open dot indicates the time convolution [cf. 
Eq. (5.9)] and the time argument of the functions A and 
B is the first (second) time argument of the Green's func- 
tion on their right (left), e.g.: 

BgBog = j dtsB{tx)g'^{H,h)B{h)g''{h,t2) 

Substituting the matrix Green's function of the form 
(5.23), we find that the collision integral becomes: 



St St 

,St^St^, 



(D3) 



(D2) 



The explicit expressions for the retarded, advanced and 
'Z' components are: 

St^ = g^(n) [B- A] + g^{n)oB{n, n)g^{n)B{n, n) 

St^ = [A;B]g'^{n) - B{n,h)g'^ {h)B{h,n)og^ (n) 

St^ = 2S{ti - t2) [A; B] - 2B{n, n)g^{n)B{h, n) 

where the (equal time) commutator is 

[A; B] = A{n, n)B{h, n) - B{n, h)A{n, n) 

The calculations performed so far are exact. However 
at one loop we are interested in terms up to the second 
order in the fluctuating fields. Then the expansion of the 
exponentials in Eq. (Dl) shows that any product of two 
functions B is proportional to terms of the form K- K- , 
which vanish after averaging over the fluctuating flelds; 
accordingly, we drop such terms. The remaining terms 
are all commutators, whose explicit (approximate) form 
is: 

[A- B] = 2i (^K_ (n) - K_ (n)) (D4) 

Note that the second order terms cancel each other ex- 
actly. The surviving first order terms lead to Eq. (5.24) 
for K_. 

For reference, we present the expression for St from 
which Eq. (5.73) and Eq. (5.76) are derived [with the 
exception of the last line of Eq. (5.76b), which is a con- 
sequence of the requirement (5.34) for K+]: 

St^=26{ti - ta) [A; B] + 2A{n, n),g^(n).4(n, n) 

-g^{n)A{n, n)A(n, n) - A{n, n)A{n, n)g^{n) 
+g^{n)B{n, n)B{h, n) + B{n, n)B{n, n)g^{n) 
+g^{n)oB{n, n)g^ {h)A{h,n) 
-A{n,h)g^{h)B{n,n)og^{n) (D5) 



APPENDIX E: DERIVATION OF EQS. (5.81b) 
AND (5.82b) 

We start the derivation by separating the Keldysh and 

^ in 

retarded/ advanced propagators contributions to Stj^ in 
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Eq. (5.75b): 

-if i fdn2dn3 3 c, n 



ley {n^r 

[2JC^(i, St, na, na) - IC^ (t + St/2, 0, ng, nz) 
-/C^(/:- 5t/2,0,n3,n2) - (n2 ^ m)] 



(El) 



Sti 
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dn2dn3 o 

712 



(E2) 



[/C^(t3,i2,n3,'n-2) - IC"{ti,t3,n2,n3) - (n2 ^ "i)] 
X [g{ti,t3, ni)g{t3,t2,n2) + 9{ti,t3, n2)g{t3, t2, ni)] 

For our convenience, we rewrote the Keldysh part in 
terms of the new time variables i, St: 



t = ^ ^ ^ , St = tl —t2, 

g{tut2) ^ 9{t,St), /C^(ti,ti) ^ /C^(f+^V2,0) 



(E3) 



Let us consider the Hmit t2 — > ii of Eq. (El); clearly 
when St ^ 0, the square bracket vanishes. However we 
know that in this Hmit g — > —2i/nSt [of. Eq. (5.31)]: 
in principle there could be a non-vanishing contribution 
from the first order expansion of the propagators in St. 
The last two Keldysh propagators depend on St in their 
first variable, but with opposite signs, and so their re- 
spective first order terms cancel each other. As for the 
first propagator, the property JC^{1, 2) = JC^{2, 1) trans- 
lates into K,^ {t, St) — IC^{t, —St), which ensures the ab- 
sence of first order terms. We conclude that in the limit 
St ^ the Keldysh propagator terms vanish. Similarly 
from the property /C^(l,2) = /C^(2,l) it follows that 
Eq. (E2) vanish for t2 = ii; this concludes the proof of 
Eq. (5.81b). 

Let us turn to Eq. (5.82b). Since dt^ — 9*2 = 2dst, 
we must expand the Keldysh propagators to the second 
order in St. At this order the square bracket in Eq. (El) 
is (up to the proper combination of angular variables): 



St^dl--^dl]lC^{t,0) 



-St"" lim dt,dt^K'^{tut2) 

*2— >tl 



where in the last equality we restored the original time 
variables. In the operator notation this is: 



se 



dtK^'dt 



Therefore 



lim (9t, -at,)Stf(ti,t2) = ^ / 
||^9t/C^9t {ti,ti,n3,ni) 



dn2dn3 3 



(ti,ti,n3,n2)| 



which proves the first part of Eq. (5.82b). 

As for Eq. (E2), using again the analytic property 

/C^(l,2) = /C^(2, 1) we conclude that when the deriva- 
tives act on the distribution functions g, the 



terms in the second line cancel each other. However there 
are non- vanishing contributions when a derivative acts on 
the propagators, such as: 



dt. 



dn2dn3 3 



lf2 9{ti,t3,ni)g{t3,ti,n2) 



{ti,t3,ni,n3) - dtK,^ (ti,i3,'«2,n3)| 



Collecting all the terms we arrive at: 

^lmi^(at, - dt^) §t^^{ti,t2) = ^ jdt3 j 



dn3dn2 q 
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(ii,i3,n2,n3) - dtK.^ (ti, i3, ^ii, ^13)} 
X [g{ti,t3,ni)g{t3,ti,n2) + g{ti,t3,n2)g{t3,ti,ni) 
which concludes the derivation of Eq. (5.82b). 



APPENDIX F: ELASTIC KERNELS IN TERMS 
OF THE INTERACTION PROPAGATOR D 

To compare the kernels in Eq. (6.9) to the correspond- 
ing expressions in Ref. 17, we use the Fourier transforms 
of Eqs. (5.51) and (5.55a) to obtain 



Re 



= —ULolm 



(Fl) 



If we assume, as is done in Ref. 17, that the Fermi liquid 
parameters do not depend on the momentum direction, 
then the interaction propagators D^'^ also do not de- 
pend on it and the above equation becomes 



Re 



-z/wim 



(F2) 



where we generalized the angular integral notation so 
that 

{C^ = J ^£^(ni,n2), n = J ^C^{n,,n2) 

We remind that our ghost propagator £9 coincides with 
the diffuson propagator D of Ref. 17. 

By substituting Eq. (F2), we rewrite the kernels (6.9) 
as: 

Sli{u;)=^S,.J^({C^){C^) - {C^n)D^ (F3) 



2d r d'^q 
V J (2^ 



(F4) 



d 



Im 



d'^q 



TTT J {2tt)'^ 



D' 



(F5) 
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The first square bracket in the kernel £ can be expressed 

as 



and using the identity: 



£3 



n; Str 



(F6) 
(F7) 



we recast it as follows: 

'5. 



(F8) 



In the second square bracket we use the identity 



vf 



(F9) 



to obtain: 



^ 

Vf 



(FIO) 



Finally the identity 

£9£9 = -idu^C^ (FIX) 
enables us to conclude that the sum of the three kernels 

that determines the correction to the conductivity - 

Eq. (6.8) coincides with the combination {Kq — Ki — 
Lq/vft) in the expression for the conductivity of Ref. 17. 



APPENDIX G: THE INELASTIC KERNEL FOR 
THE PHASE RELAXATION TIME 

Let us consider a uniform system in which the bosons 
are assumed to be in equilibrium with the electrons. In 
other words, the distribution function / is independent of 
r, n and the boson-electron collision integral (3.18) must 
vanish. The latter condition enables us to express the 
bosonic distributions N'^ in terms of / and to obtain 

Tij(e,w) = -QdSiuiTij)^ j dei'^{e,ei-,u!) (Gl) 

[from now on irrelevant angular and momentum variables 
are omitted; all relevant definitions can be found in Sec. 
Ill] . The former condition implies that the collision inte- 
grals (3.19d). (3.20b) and (3.20f) vanish and therefore the 
kinetic equation for (the zeroth harmonic of) / reduces 
to 



dtf{e;t)= dw deiA{uj)^{e,ei;u;) (G2) 



where 

A(co) =^i^ReTr( hstrt^ + 1 



St^ 



(03) 



Wc substitute Eq. (Fl) and the similar relation for the 
triplet channel [D^ being the triplet channel propagator] 
into the expression for A{co); then we use the identity 
(5.44) and obtain 



A{u)) = — ^ImTr \strb (d^ + dJ^) b 



(G4) 



Using again Eq. (5.44), we recover immediately the form 
of the inelastic kernel Aioj) given in Ref. 27. 
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formed independetly by M.Yu. Reyzer (private communi- 
cation). 

It was shown by G.V. Chester and A. Thellung, that 
Eq. (1.1) remains valid for arbitrary scattering strength 
as long as the scattering rates and the density of states are 
smooth (C2) function of energy near the Fermi level. 
^® The quantum kinetic equation with the necessary conser- 
vation laws was not actually derived in Ref. 20, so we will 
not be able to compare their approach with ours. 



^"^ As given in Eq. (3.1), the equilibrium distribution func- 
tions N"q{uj) are defined only for a; > 0; for uo < Q they are 
found by using the property in Eq. (3.2). 
The current operator used in Ref. 21 does not satisfy the 
continuity equation for the long range interaction poten- 
tial. 

^® These solution are exact only in the absence of the elec- 
tric field, since in its gauge invariant form the operator 
acting on Sg± is the same that appears on the left hand 
side of Eq. (5.71). We could include perturbatively field 
dependent corrections to our expressions, which would be 
of the first order in E for 5g- and of the second order for 
5g+ [since 5g+ vanishes in equilibrium] . However, as noted 
before, the first property in Eq. (5.79) implies that these 
corrections cannot contribute to the physical quantities we 
are interested in and therefore we do not include them in 
our calculations. 

""^ The operator M is gauge-invaxiant and it has the same 
form in terms of g or g^ . 

'^^ by the simple substitution C - C ^ C + - AC in 
the kernels. 

The Boltzmann equation description of strictly one- 
dimensional systems is not applicable and considering the 
quasi one-dimensional ballistic case within our scheme is 
meaningless, because of the effects of boundary scattering. 
The ballistic regime in three dimensions also can not be 
considered within our scheme, because the main effect on 
the thermal conductivity is due to the inelastic scattering 
processes with momentum transfer of the order oi kp- 
We choose to collect the bosonic contributions into a single 
kernel so that the resulting momentum integral is conver- 
gent. 

They can be obtained from Eqs. (6.25)-(6.26) in the limit 
d ^ 2. 

A similar argument is presented in Ref. 21. 
This definition of the density of states is one half of the one 
given in Eqs. (2.19), because of the different limits for the 
oj integration in the energy density and the specific heat. 



